Incompressible Limit of the Compressible Hydrodynamic Flow 

of Liquid Crystals 

o 

u . 

Q*. Shijin Ding* Jinrui Huang''" Huanyao Wen* Ruizhao Zi§ 

<N 

i 1 Abstract 

Ph 

This paper is concerned with the incompressible limit of the compressible hydro- 
, dynamic flow of liquid crystals with periodic boundary conditions in WL N (JV = 2,3). 

It is rigorously shown that the local (and global) strong solution of the compressible 
system converges to the local (and global) strong solution of the incompressible system. 
Furthermore, the convergence rates are also obtained in some sense. 
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1 Introduction 

The hydrodynamic flow of incompressible liquid crystals was first derived by Ericksen [6] 
and Leslie [16] in 1960s. However, its rigorous mathematical analysis did not take place 
until 1990s, when Lin |17j and Lin and Liu [18 } 119 } 120] made some very important progress 
towards the existence of global weak solutions and partial regularity of the incompressible 
hydrodynamic flow equation of liquid crystals. 

In the context of hydrodynamics, the basic variable is the flow map (particle trajectory) 
x(X, t). X is the original labeling (Lagrangian coordinate) of the particle. It is also referred 
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to as material coordinate, x is the current (Eulerian) coordinate and referred to as reference 
coordinate. For a given velocity field u(x,t), the flow map is defined by the following 
ordinary differential equation: 

x t (X,t) = u(x(X,t),t), x(X,0)=X. (1.1) 

The deformation tensor F(X, t) is defined as 

(1-2) 

Applying the chain rule, we see that F(x, t) satisfies the following transport equation: 

,_, d „ dxt du du dx _ _, ,„ „. 

F t + u • VF = — F = —± = — = — - — = VuF, (1.3) 
dt dX dX dx dX ' v ; 

which stands for F^ + UkVkFij = VkViF^j- 
Let us define the density as 

(1-4) 

By the identity of the variation of the determinant of a tensor 

JdetF = det Ftr(F~ 1 5F), (1.5) 

we have 

then we get the transport equation 

d t p + u- Vp + P V -u = 0, (1.7) 

which can also be reached by the conservation of mass. 

The momentum equation of motion for hydrodynamics flow of liquid crystals can be 
derived from the least action principle (Hamilton principle). The action functional takes 
the form: 

A(x) = jT Q po (x)[x t [ 2 - W x {F)^jdXdt, A>0, (1.8) 

where is the original domain occupied by the material. The first part of A{x) denotes 
the kinetic energy and the second one denotes the elastic energy. By the definition of the 
density in (|1.4p . we get 



p (X)\x t \ 2 dXdt = [ [ ^^\x t \ 2 det FdXdt = [ [ p(x, t)\u(x, t)\ 2 dxdt. (1.9) 



o Jn Jo Jn detF j j n 



In this paper, we will consider the isotropic energy functions W (F) of the form 

W\F) = ( X 2 uj ( + V - \F-'V x n (X)\ 2 ) detF, (1.10) 



^ det F J 2 

where u is a C°° function and denotes the energy density, uq(X) : — > S 2 is a unit-vector 
field that represents the molecular orientation of the liquid crystal material. The first term 
of the R.H.S. of (jl.lOp should be regarded as a penalization term which drives the motion 
toward incompressibility in the limit as the parameter A becomes large. We can use the 
kinematic assumption to give 

n t +u-Vn = 0, (1.11) 

then we have 

n(x(X,t),t)=n (X), (1.12) 

which implies that the center of gravity of the molecules moves along the streamline of 
the velocity, and the second term of the R.H.S. of (jl.lOj) | \F^ 1 Vxno(X)\ 2 is equivalent 
to ||Vn| 2 in Eulerian coordinate, where v > is viscosity of the fluid and denotes the 
microscopic elastic relaxation time. 

For any y = y 2 , • • • ,Vn) G C%(n-,R N ), let x £ = x + ey(x), F e = and p e = 

de le=0 



p(x £ (X, t), t). Then we apply 5 = 4- \ =n to (1.8) and get 



= 5A(x € ) =11 \ Po (X)(x t ,y t ) - 5W X (F)] dXdt, (1.13) 
Jo Jq 1 J 

where we denote the inner product / • g by (f,g) for some /, g : O — > R . By using the 

definition of p(x,t), we have 

/ p (X)(x t ,y t )dXdt = - If po(X)(x tt ,y)dXdt = - f [ ^l{ Xtu y)dxdt 
o Jn jo Jn Jo Jn aet±l 

= - I I p(x,t)(u t + u -Vu,y)dxdt. (1-14) 
Jo Jn 

Now we turn to show the calculation about the second term of the R.H.S of (|1 . 13|) . 
Noting the facts that 



det F \ dx dX ) det F 

and 



(5 (det F £ ) = det Fir (F~ x V x y) = det FV • y, (1.16) 



we use the definition of p to give 



-A 2 f I 5 (uf I an r 1 dXdl 



o ^n 



det F e 

-A 2 f / f-^'f^^^detF + c.f^ldetFV-yW/AA// 



Jflo 



det F / det F V det F 



-A 2 / / (uj(p) - poj'(p))V -ydxdt 
Jo Jn 

-A 2 / / (VP(p),y)dxdt, (1.17) 



where we define the pressure as 

P(p) = -o;(p) + pa; / (p) ! (1.18) 

which can also be reached by the first law of thermodynamics. Here we show the proof 
briefly. According to the first law of thermodynamics, we obtain 

dW = -PdV - SdT, (1.19) 

where W, V, S and T denote the energy, volume, entropy, and temperature, respectively. 
Moreover, 

in 

W = loV, V=—, oj = uj(p), (1.20) 
P 

where m is the total mass. From (|1.19p . we have 

„ dW d(uV) Tr dto r r du} (w) ,m , 

then we get (1.18). 

On the other hand, by the Leibniz's formula, we have 

--I t S (UF^V x n (X)\ 2 det F e ) dXdt 
2 Jo Jn ^ ' 

- V - f [ [2(F- 1 V x n (X),d(F^- 1 V x n (X)) det F + \F~ l V x n (X)\ 2 5 (det F e )] dXdt, 
1 Jo Jn 



= V 1 + V 2 , (1.22) 
where 

V 2 = ~ I [ \F- Y V x n {X)\ 2 det Ftr (F _1 £F e ) dX dt 
2 Jo Jn 
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-ff \Vn\ 2 detFV ■y—^—dxdt 

- / / |Vn| 2 V • ydxdt 
2 Jo Jn 

[ (V(\Vn\ 2 ),y)dxdt, 
o Jn 



(1.23) 



and 



I f (F- 1 Vxn {X),S(F e )- 1 Vxn {X)) det FdXdt 
Jo Jn 



o Jn. 



dX 

d Xj n (X)-^-,d Xr n (X)S 



dX r 



det FdXdt 



dX r d dx% dXi 
dxi de dXi dx e - 



~ v t L ( d ^M x )^-' d x r Mx) 



dXdt 



= V 



o Jn 

t 



dy k \ 



1 ' k dxij detF 



dxdt 



v\ f (V • (V 

jo Jn 



n <S> Vn) , y) dxdt, 



where we differentiate F e (F € ) 1 = Ijv with respect to e and obtain 

'dX r 



Ox" 



dX r d dx% dXi 
dx e k de dXi dx\ 



(1.24) 



(1.25) 



t=0 



Particularly, for the case of the incompressible system, the calculation above can be simpli- 
fied by the incompressibility, i.e. detF = 1. 

Finally, we turn to deal with the dissipation. We consider the case 

d 



dt 



j^total ^ 



(1.26) 



where E total denotes the summation of kinetic energy and the elastic energy shown as follows 



E 



total 



p\u\ 2 + \ 2 S{p) + -\Vn\ 2 )dx, 



and A denotes dissipation as follows 

A = / (2fi\D\ 2 + k\V ■ u\ 2 + 6\n t + u ■ Vn\ 2 ) dx, 
Jn 



(1.27) 



(1.28) 



wh ere D = v "+ vT " = f v ^+ v ^ 



NxN 



is so-called deformation tensor, \x and k are the 



shear viscosity and the bulk viscosity coefficients satisfying p, > and Nk + 2fi > 0. Let 



u T = u + tv, v e Clip,; R N ). Then applying 6 = 
parts, we have 



to A, and by using integration by 



T=0 



SA 



p 



(Vu + V T u, Vv + V T v) dx + K / (V ■ u, V ■ v)d. 



= 77 / [(V«,Vu) + (V T u,V T u)] dx + ^ / [(Vu,V T u) + (V T u,Vv)] da; 
2 Jo ^ Jn 

+k / (V • it, V ■ t>)d:r 
Jo 

= / (pAu + (k + m)V(V •«),«) ds = 0, (1.29) 
Jo 

where the third term of ()1.28[) denotes the microscopic dissipation, and has been treated 
as zero when we apply 5 to (jl.28p . Clearly, if the system is incompressible, the term 
(k + /u)V(V • u) in (fL29l) should be zero. 

Combining (j!.14j) . (|1.17p . (jl.23[) . (|1.24j) and (jl.29p . we have the momentum equation 

/ |Vnl 2 \ 

put + p(u ■ V)u + A 2 V (P(p)) = pAu + (« + /u)V(V • it) - uV ■ Vn <g> Vn - - ^2jv .(1.30) 



2 

Finally, we deal with the heat flow for harmonic maps. For any C 1 map n : £1 — > S 2 , 
6 f 

let I(n) = — J \\7n\ 2 dx denote the energy of n, where 9 > is the viscosity constant. The 
2 Jo 

critical points of the energy I(n) are called harmonic maps. To get the critical points of 
the energy I(n), let 

n + rip 

mm = i r- (1-31) 

' \n + Tip\ y ' 



Then by applying 5 = 4- 



to J(m(r)), we have 



T=0 



= -5 ( I \Vm( T )\ 2 dx J = / Vm(0) : Vm'(0)dx = / Vn : Vm'(0)dz. (1.32) 
2 VJo / Jo Jo 

Direct calculation from (jl.3ip yields 

m'(r) = , ^ , - (n + ^ )(n+ ,7^ ) , (1.33) 
|n + Tip\ \n + ry| 

and by using the fact that \n\ = 1, we obtain m'(0) = <p — (n,ip)n. Consequently, we have 

/ Vn : Vm'(0)dx = / Vn : V[y — (n, y?)n]dx = / Vn : \7ipdx — / Vn : V[(n, </?)n]da; 
Jo Jo Jo Jo 

(An,(^)dx- / Vn : V[(n, v?)n]dx. (1.34) 



Noting that V[(n, </?)n] = (n, ip)Vn + (V<^ • n + Vn • y?)n, one can easily deduces that 

Vn : V[(n, </>)rt] = |Vn| 2 (n, tp), (1.35) 

where we have used the fact V (|n| 2 ) = 0. 

Combining (TQ21 . (TOP and (fOoD . we have 

/ (0(An+ |Vn| 2 n),v9) = 0, (1.36) 
Jo 



then we get the heat flow for harmonic maps by considering the evolution equation associ- 
ated with the harmonic maps 

n t + (it • V)n = 6>(An+ |Vn| 2 n). (1.37) 

In conclusion, the compressible hydrodynamic flow equation of liquid crystals in = 
can be written as follows (the functions and the viscosity constants should depend 
on the value of the parameter A): 

p\ + V • (p A it A ) = 0, 

(p x u x ) t + V • (p x u x it A ) + A 2 V(P(p A )) 

(1.38) 

= p x Au x + {k x + p x )V(V ■ u x ) - u x V ■ (Vn x Vn A - ^p^I N ) , 
n A + (u x ■ V)n A = 9 x (An x + |Vn A | 2 n A ). 

Here x € T^, a torus in R N , N =2 or 3, t > 0, the coefficients p x ,v x ,6 x > 0, and 
2p x + Nk x > 0. The unknown functions are the density p x , the velocity u x € R , and the 
molecular orientation of the liquid crystal material n x : T N x (0, +oo) — > S 2 , which is a 
unit vector, i.e., |re A | = 1. P(p) is the smooth pressure-density function with P'(p) > for 
p > 0. And the symbol <g> and denote the tensor product such that u®u = (u l u J )i<i t j<N 
and Vn Vn=(n Xi ■ n Xj )i<i t j<N, respectively. 

From the mathematical point of view, it is reasonable to expect that, as p x — > 1, the 
first equation in (jl.38p yields the incompressible condition V • it = 0. Suppose that the 
limits u x — > u and n x — > n exist as A — > oo, and the viscosity coefficients satisfy that 

p, x ->• p. > 0, K X — y k, z/ -> v > 0, 6> A -)• 6» > 0, as A -> oo, (1.39) 

then (at least formally) we obtain the following incompressible model by taking A — > oo: 

V • it = 0, 

' ut + (u ■ V)u + Vp = /uAu - t^V • (Vn Vn), (1-40) 
nt + (u • V)n = 6(An + |Vn| 2 n), 

where p is the limit of the term prV (P(/3 A )) — ^ A ^ V ^ ' • For simplicity we assume in the 
following that k x = k, p x = p, p x = p and X = 6 are constants independent of A, satisfying 
that p, v, 6 > and 2p + Nk > 0. 

In a series of papers, Lin [T7] and Lin and Liu [181 EE I2D] addressed the existence and 
partial regularity theory of suitable weak solutions to the incompressible hydrodynamic 
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flow of liquid crystals of variable length. More precisely, they considered the approximate 
equation of incompressible hydrodynamic flow of liquid crystals (i.e., |Vn| 2 is replaced by 

(]_ — 1 77, 1 ^ 77 

2 ), and proved [18) . among many other results, the local existence of classical 

solutions and the global existence of weak solutions in dimension two and three. For any 
fixed e > 0, they also showed the existence and uniqueness of global classical solution 
either in dimension two or dimension three when the fluid viscosity [i is sufficiently large; 
in [20], Lin and Liu extended the classical theorem by Caffarelli-Kohn-Nirenberg [1] on the 
Navier-Stokes equation that asserts the one dimensional parabolic Hausdorff measure of the 
singular set of any suitable weak solution is zero. See also [22} 1231 [28] for relevant results. 

The paper is organized as follows. In Section 2, we state the main results of this paper. 
In Section 3, we prove that the strong local solutions of (jl.38p exist for sufficiently small 
disturbances from the general incompressible initial data, meanwhile, we get the uniform 
stability of the local solutions family which yields a lifespan of the system (|1.38j) . Based 
on it, we show that the local solutions of (jl.38p converge to a local solution of the limiting 
incompressible system (jl .401) by means of compactness arguments. In Section 4, the global 
existence of the strong solution to the incompressible system (ll.40p is drived provided the 
initial data of the incompressible model sufficiently small. In Section 5, we obtain the 
convergence rates about (p x ,u x ,n x ) — > (l,u,n) in some sense when A — > oo. These results 
depend on some methods modified from [11] \13\ [T5l 127] . 

2 Notations and Statements of Main Results 

Noting that 

/ I Vr? A l 2 \ 

V ■ (Vn x QVn x - ] -^^I N ) =Y J An i^n x , (2.1) 

we can rewrite (11 .38H as follows when \p x — 1| is small: 
Pt+V- (p x u x ) = 0, 

< u x + (u x • V)u A + ^VP(p x ) = f x Au x + ^V(V • u x ) - f Eili An A Vn A , (2-2) 
n x + (u x ■ V)n x = 9(An x + |Vn A | 2 n A ). 

Throughout this paper, for simplicity, we will denote Y2i=i An x Vn x by An x ■ Vn A , the 
norms m L 2 (T n ), H s (T n ) and L°°(T n ) by || • ||, || • || s and || • respectively, and especially 
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the norms in L (T n ) by || • ||^4. We define 

E.(U(t)) = hT, [ [\ 2 \V a (p-l)\ 2 + \V a u\ 2 + \V a n\ 2 ]dx, 

\a\<sJT n 



Es{U{t)) = \ E 



x\< S JT 



A" 



-\V a (p - 1)|" + p\V a u\ 2 + |V a n' 2 



P 



(2.3) 



where U = (p,u,n). We point out that the parameter A is hidden in the definition of 
E s (U(t)) and E s (U(t)). Similarly, we define 

W))4E / [X 2 \V a (p-l)\ 2 + \V a u\ 2 ]dx + l V / |VV^n| 2 dx, 



l/3|<s-l' 



^(^W) = I E 



A 2^P)| va(p _ 1) |2 +p | V a n j2 



dx + - V / |W^| 2 dx, 



(2.4) 



|/5|<s— 1 



where there are some differences between (12.31) and (12.41) in terms of n. It is obvious that 



(2.5) 



E s (U(t)) ~ £7,(I7(t)), F s (£/(i)) - F s (C/(t)), 



provided that |p — 1| is sufficiently small. 

Now we state the main results of this paper. 



Theorem 2.1 Consider the compressible model il.38\) with the following initial data 

n (x) + ng(x) 



p A (x,0) = l + p$(x), u A (x,0) = u (x) + u A (x), n A (x,0) 



\n (x) + n$(x)\ ' 



(2.6) 



where Uq, no satisfy 

u (x) e H s+1 (T n ), V-u = 0, n e H s+1 (T n ,S 2 ), (2.7) 



/or any s > [4j-] + 2. Moreover, for some small positive constants 5o, the functions Pq(x), 
Uq(x), Uq(x) are assumed to satisfy 



i— A 



< A- 2 5 , ||u A (x)|| s+1 < A-^o, ||Vn A (x,0) - Vn (x)|| s < A" 1 ^. 



(2.8) 



Then the following statements hold: 

Uniform stability: There exist fixed constants Tq and C independent of A such that the 
unique strong solution (/3 A ,n A ,n A ) of system A1.38\) exists for all large A on the time interval 
[0,Tq]. Furthermore, the solution family satisfy: 



'o 



A I 2 



<c, 



t r 



E s ^(d t U\t)) + / MllV^ll^i + ^ + ^IIV-^H^ + eilV^ll^i <C, (2.9) 
Jo L J 



|n A | = 1, in Q Tq = T n x [0,T ], 



where E s _ l (d t U(t)) = \ £ f Jn (X 2 \V^d t p\ 2 + \V^d t u\ 2 + \V^d t n\ 2 ) dx. 

W<s-1 

Local existence of solutions for incompressible system: There exist functions u and n 



such that 
( 



p x ^l in L°°{[0,T o Y,H')nLip{[0 ) T o Y,H'- 1 ) t 

(u x ,n x ) -» (u,n) weakly* in L°°([0,T ]; H s ) D Lip([0,T ]; H 3 - 1 ) 

(u x ,n x ) -> (u,n) in C([0,T }; H s '), 



(2.10) 



for any s' in [0,s), and the function pair (u,n) is the unique strong solution of the incom- 
pressible system of liquid crystal 1 1.40\ ) with the initial data 



u(x,0)=uq, n(x,Q)=riQ, 



(2.11) 



for somepe L°°([0, T \; H 8 ' 1 ) n L 2 ([0, T \; H s ). 

Theorem 2.2 Consider the strong solutions (p x ,u x ,n x ) of system (1.38\) obtained in The- 
orem \2J\ Suppose in addition that the initial data satisfies 



Nils + l|Vn ||s-i < eo, 



(2.12) 



where Eq is a positive constant. If Eq and A -1 are sufficiently small, then for every fixed 
T > 0, the strong solution (p x ,u x ,n x ) exists on the time interval [0, T x ) with T x > T, which 
satisfies 

F s {U x (t)) + / [p\\Vu\\ 2 s + (k + p)\\V ■ u\\ 2 s + 0||V 2 n||2_i] < 4(e + \~ 2 Sq), (2.13) 
J 

and 



E s ^{d t U x {t))+ / mIIV^II^i + Ck + ^IIV-^H^! +6\\Vdtn A \\i-i < CexpCt, (2.14) 
Jo 1 J 

for t G [0,T], and T x — > 00 as A — >• 00. Furthermore, as A — >■ 00, (p x ,u x ,n x ) converges 

to the unique global strong solution (l,u,n) of the incompressible system of liquid crystal 

{l.^Cfy , with \n{x,t)\ = 1 in Qt and 



Mis + l|Vn||?_i + / (p\\Vu\\ 2 s + ff||V 2 n||2_i) < C* l£o , 



(2.15) 

where C\ is a uniform constant depending only on N, s, p, k, v, 6 and the domain, but 
independent of £q and A. 



m 



Theorem 2.3 Under the assumptions of Theorem \2.1l the convergence rate of p x , u x and 
n x when A — > oo is deduced in the following sense that 

X\\p x - 1|| 2 + \\u x - u\\ 2 + \\n x - n\\ 2 2 + J [\\u x - u\\\ + \\n x - n\\fj < CX' 1 . (2.16) 

Furthermore, we have 

||V(/-1)||2_ 2 <C7A- 4 . (2.17) 

The statement also holds for the global strong solution given in Theorem \2.2[ 

Remark 2.1 In Theorem \2. 31 we do not give the convergence rates about the higher order 
derivatives of u x and n x because we know little about the convergence rate of the pressure. 

3 Local Existence and Uniform Stability 

Let Uq = (l+/jg , no+Ug , i" ^"" i )• We consider a set of functions Bj, (Uo) = Bh, (Uo, 5, s, K\ , K2) 

contained in L°°([0, T ];H S ) n Lip([0, T ]; IP -1 ) with s > 3 and defined by 
r 

\\{jp - 1)| + \u - u \ + |Vn - Vn | < 5, 
< E s (U(t)) + / [m||V«||2 + (k + /x)||V • + 0||Vn||2] < (3.1) 

t 

E s ^(d t U(t)) + / [/i||V^tu||2_i + (« + M)||V • a^H^ + 0||V$n||2_i] < if 2 . 
l> JO 

For V = (£ A ,i> A ,m A ) G ^ o ([/ ), define U = (p x ,u x ,n x ) = A(V) as the unique solution of 
the following problem: 

p\ + (v A • V)p A + £ A V • u x = 0, 
< u A + (w A • V)u A + A 2 ^|^Vp A = j^Au x + to) V(V • u A ) - ^An A • Vn A , (3.2) 
n\ + (t; A • V)n A = 6>(An A + |Vm A | 2 n A ). 

We plan to show that for appropriate choices of To, 5, K\ and K2 independent of A, A 
maps Bj, (Uo) into itself and it is a contraction in certain function spaces. We emphasize 
that the solutions will depend on the value of the parameter A, but for convenience, the 
dependence will not always be displayed in this section. 

We need the following lemma for the proofs. 

Lemma 3.1 \13\ Assume that f,g£ H"(T ). Then for any multi-index a = (a±, . . . ,ayv)> 
oti G N, [a| < k, we have 

(i) \\D a (fg)\\ < C N>k (\\f\\oo\\D k g\\ + \\g\U\D k f\\), 
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(ii) \\D a (fg) - fD»g\\ < C N , k (\\Df\U\D k - l g\\ + \\g\U\D* f\\) 



where C^k is a constant depending only on k and N. 



Before proceeding any further, we apply D a to (|3.2|) and get 
d t D a p + (v ■ V)D a p + £V • D a u = ni, 

d t D a u + (y ■ V)D a u + X 2 ^VD a p = |AD Q u + ^V(D Q V • u) 
d t D a n + D a (v Vn) = 8AD a n + 9D a (\ Vm| 2 n), 

where 



|D Q (An- Vn) + n 2 , 



ni 
n 2 



-[£> Q (t> • Vp) - v ■ VD a p] - [ZT(fV ■ u) - £V ■ 
-[D Q (t> • Vn) - v ■ VD a u] - A 2 



+ 



D Q ( ^Am) - —AD a u 



V(V • «) 



K + p 



VD a (X7 -u) 



D' 



7 (An • Vn) 



-£> a (An • Vn) 



We will prove that A maps B^, (Z7o) into itself by two steps and denote by C the constants 
independent of A, K\ and K 2 in these two steps. 
Step One: Estimates about n. 

Taking the L 2 inner product of (|3.3p ^ with D a n, and using integration by parts, we have 



1 £ I \D a n\ 2 + 6 [ \VD a n\ 2 

2 dt Jjn Jfn 

-[ D a {v -Vn) ■ D a n + 9 [ D a (\Vm\ 2 n) ■ D a n 



N 1 +N 2 . 



(3.4) 



Next, we will give the estimates about N\ and N 2 in the following three cases respectively. 
Case 1: \a\ = or D a = d t . 

When \a\ = 0, we use integration by parts and the Sobolev embedding H 2 (T N ) 
L 00 ^) for N = 2,3 to get 



\Ni\ 



(V-«)^- 



< ^l|V-v||oo||n|| 2 < CK\\\n\\ 2 



(3.5) 



and 



|JV 2 | = I |Vm| 2 |n| 2 cfo; < 0||Vm|£,||n|| 2 < CK x \\nf. 
Jt n 



(3.6) 
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On the other hand, when D a = dt, we use the Cauchy inequality to give 



|iVi 



/ (vf Vn) -n t - (v ■ Vn t ) 
Jt n Jt n 



< IKIIoollVniiKH + iHiooiiVnilllKil 

< gl|Vn t f + C7(||Vn|| 2 + H|l)||n t || 2 + C|HH 

< ?l|Vn t || 2 + C(\\Vnf + K^Wntf + CK 2 



(3.7) 



and 



20 I (Vm:Vm t )(n-n t ) + e \Vm\ 2 \n t \" 

fN JfN 

2 ll„ ||2 



l^| = 

< 2#||Vm|| 00 ||Vm t ||||n|| 00 ||n t || +9\\Vm 

< C (||Vm|| 2 ||Vm t || 2 + ||Vm|||) ||n t || 2 + C\\n\\ 2 2 

< CiK^ + K 1 )\\n t \\ 2 + C\\n\\l 

Case 2: D a = ViV? with \/3\ < s - 1. 

Integrating by parts and using Lemma 13.11 we have 



(3.8) 



liVi 



/ V p (v ■ Vn) 
Jt n 



n 



< C\\V p {v- Vn)||||AV^ra|| 

< Cdl^HoollVnll + HVnlUlV^IDIlV'V^nl 



< 7l|V 2 V /3 n|| + CK4V 



s-l- 



(3.9) 



|V 2 V /3 



n 



Similarly, we have 

|jV 2 | = -9 [ Vm| 2 n) • AV^n 
Jt n 

< C||V^(|Vm| 2 n)||||AV^n|| 

< C (HVmll^UV^nll + UnlUHV^dVml 2 ) 

< C (HVmll^llV-^H + llnlloollVmlloollVmll) ||V 2 V^| 

< e -\\V 2 vPn\\ 2 + CK 2 \\n\\ 2 s . 

Case 3: D a = VjV^ with | 7 | < s - 2. 

It follows from Lemma 13.11 and the Cauchy inequality that 



(3.10) 



< c 



[ vP(v t ■ Vn) • V p n t - [ V p {vVn t )-V p 
Jt n Jt n 

\\V^v t - Vn)|| + HV'VVnt)! 



N p n t \ 



< C7 (||v*||oo||V s ra|| + HVnlUIV^VH + |MU|V s n t || + HVratlloollV'-^H) |M| s -i 
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< ^rH\\ 2 a - 1 + C(K 1 +K 2 )\\n t \\ 2 s - 1 + C\\n\\l (3.11) 
On the other hand, noting that 

N 2 = -20 V 7 [(Vm : Vm t )n] • AV 7 n 4 - / V 7 (|Vm| 2 n t ) • AV 7 ra t 

JjN JfN 

= N 2A + N 2<2 , (3.12) 
and using Lemma 13.11 and the Cauchy inequality, we have 
\N 2)2 \ < C||V 7 (|Vm| 2 n t )||||AV 7 m,|| 

< adiVmii^nv-^ii + Kiloollv-^ivm^iDiivntiU-! 

< C (|| Vm||i||V s - 2 n t || + ||n t || 2 ||Vm|| 0O ||V s - 1 m||) \\Vn t \\ s ^ 

< gllVntll^ + C^llntll^!. (3.13) 

Nevertheless, the estimates of N 2 ,i are slightly different between the cases s = 3 and s > 4. 
Indeed, if s > 4, then 

[AT 2;1 | < C||V 7 [(Vm : Vm t )n]||||AV 7 n t || 

< C (llVmlloollVmtHoollV^nll + ||n||oo||V a - 2 (Vm : Vm t )||) \\\\Vn t \\ s ^ 

< C (||Vm|| 2 ||Vm t || 2 ||V s - 2 n|| + ||n|| 2 || Vm\\ 0O \\V'- 1 m t \\ 
+ ||n|| 2 ||Vm t || 00 ||V s - 1 m||) ||||Vn t || s _i 

< ^WMl-l + CK^Wnf^, (3.14) 

and if s = 3 (D a = Vjdt or VjVjSt), then we use another Sobolev embedding H 1 (T N ) 
L*(T N ) to get 

|iV 2) i| < C||V 7 [(Vm : Vm t )n] \\ || AV 7 m.|| 

< C{\\ (Vm : Vm t )n\\ + \\Vj [(Vffl : Vm t )n] ||} ||Vn t || s _i 

< C [|| (Vm : Vm t ) n\\ + || (VjVm : Vm t ) n\\ 

+ || (Vm : VjVmt) n\\ + || (Vm : Vm,) Vjn\\] ||Vnt|| s -i 

< C (||Vm||i4||Vm t || i 4||n|| 00 + ||V 2 m|| I/ 4||Vmi|| 

L 4 ll^lloo 

+C||Vm|| 00 ||V 2 m t ||||n|| 00 + || Vm|| L 4 1| Vm t || Vn||oo) ||Vn t || s -i 

< C(||Vm|| 2 ||Vmi|| 1 ||n|| 2 +C'||Vm|| 2 ||V 2 mi||||n|| 2 + ||Vm|| 1 ||Vmi|| 1 ||Vn|| 00 ) ||Vn t || s _i 

< ^llVntll^i + C^x^Hnll 2 . (3.15) 



14 



In conclusion, substituting the estimates (|3.5p . (|3.6p . (|3.9p and (|3.10p into (|3.4p . and 
using the Cauchy inequality, we obtain 

|E HV Q n|| 2 + ^ ||W«n[| a <C(if? + l)H2. (3.16) 

a|<s |ct| <s 

Owing to |n(x,0)| = |rto(x)| = 1, x G T , we have ||n(x,0)|| = ||no(x)||. Together with the 
Gronwall's inequality, (|2.8p and (|3.16p . we find that 



n 



|2 < e C(^+l)T || n(a , )0) ||2 < e C(Kl+l)T Q| no ||2 + ^2^2) < ^ (3J7) 



where we have chosen Tq small enough for Tq < T\ = } and A > 1. Furthermore, 
integrating (|3.16p over [0,i] and using (|3.17p . we obtain 

f \\Vn\\l < C [Kl + 1) T + ||n || 2 + A~ 2 <5 2 < C. (3.18) 
J o 

On the other hand, substituting the estimates fl377|), §3Ji$) and (pmT) - (j3TT5l) into 
and then using (|3.17p and the Cauchy inequality, we obtain 

| E / W + e E / I 

\P\<s-l " |/3|<«-1 " 

< C^ + tff + l) |H|2_i + C(if? + ^ 2 2 + l). (3.19) 
Recalling that the constrains on the initial data and (|1.38p ^. we use Lemma 13.11 to get 

|Mx,0)|| s _! < C [\\(u(x,0) ■ V)n(s,0)|| a _i + ||An(x,0) + \Vn(x, 0)| 2 n(x, 0)|| s _i] 

< C (||«o|| a -i + A- 1 ^) (||Vn Q || s _i + A- 1 ^) + C (||Vn || s + A" 1 ^) 

+C (HV^ollS-i + A-^o 2 ) (||Vno||.- 2 + X-%) + C (|| Vn ||«-i + A^ 2 ) 

< C. (3.20) 

Then by the Gronwall's inequality , we have 

_! < Ce C K 2+ ^ +1 ) T ° [IK^O)^.! + + K 2 + 1)T ] < C, (3.21) 



2 



provided that To < Ti = ^-p- . Furthermore, integrating (I3.19P over [0,t] leads us to 

^ /"* IIVthK?-! < C. (3.22) 
J 

Step Two: Estimates about p and u. 

Taking the I? inner product of the first two equations of (|3.3p with X 2 P ^ D a (p — 1) 
and ^D a u respectively, and using integration by parts, we arrive at 



1 d 
2~dl 



^p-\\D a (p-l)\ 2 +Z\D a u\ 2 
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a u\ 2 



+fi [ \VD a u\ 2 + (n + p) [ |V--D 

Jfn Jfn 

h+h + h + h + h + h + h + h + h, 



(3.23) 



where 

h = 

h = 

h = 

h = 

h = 

h = 

h = 

h = 



1 

2 JT 

1 

2 j T 

2 



\XD a (p-l)\ 2 d t 1 p ^ * ' lrv 2 



\\D a (p - 1)| 2 V • 



+ tt\D a u 
v ] + V- (£v)\D a u\ 2 



X 



JT" 

D a (p - 1)( [D a {v ■ Vp) - v ■ VD a p] + [D a {£V • - £V • D a u] }, 



_ A 2 / P'(0 ^ 



v [ D a (An • Vn) • D a u, 
/ £[D a (v -Vu) -v -VD a u]- D a u, 

JT" 



-a- / e 

/T" 



{ 



T™ 



D a I ^Au 



+ 



/ ({ 


D a 


Jjn 








Jjn 





K + p 



AD a u 
V(V • «) 



_ K + p 



D a u, 



VD a (V ■ u) \ ■ D a u, 



-(An • Vn) 



-D a (An -Vn)}- D a u. 



We will estimate all these terms in the three cases as mentioned in Step One. 

First of all, we give the estimates of I1—I3 for all the three cases. Without losing any 
generality, we assume that A > 1. Furthermore, we choose 5 small enough so that (|2.5[) and 
|£ — 1| < I hold. Then by the smoothness of the pressure P(-) and the Sobolev embedding 
H 2 (T N ) ^ L°°(T N ), we have 



I'll < 



a..\\2 



Ut\\ 00 \\>jy'(p-i)\f + Ut\u\D a u 



< G||6|| 3 (||AD°(p-l)|| a + ||i3 a u|| a ) 

< C\~ x k\ (\\XD a (p - 1)|| 2 + \\D a u\\ 2 ) , 



\h\ < 



£p"(0 - p'(0 



p'(0 



llVfUooHoo + 

+ (||Ve||oo|b||oo + U\U\V-v\U\\D a U \\ 2 

< C (A- 1 ||A(<e - l)|| 3 ||«|| a + \\v\\ 3 )(\\\D a (p - 1)|| 2 + H^nll 2 ) 

< c(\- l K 1+ Kf] (\\XD a (p-l)f + \\D a u \\ 2 ), 



(3.24) 

V-r|U ) \\XD a (p-l)\\ 2 



(3.25) 



1fi 



\h\ < A||P"(Olloo||Ve||oo(||Ai? a (p-l)|| 2 + ||£> a «|| 2 ) 

< C||AVeib (!|A^(p - 1)|| 2 + ii^^ii 2 ) 

< CK\ (II XD a (p-l) || 2 + \\D a u\\ 2 ) . (3.26) 

Secondly, we give the estimates about -Z4-/9 when |a| = or D a = dt- When |a| = 0, 
the quantities Ij, 4 < j < 9, are all equal to zero except for 1$. Thus we need only to 
estimate J5. By using (|3.17|) and integration by parts, we have 

141 < C||«||||Vn||oo||An|| < C\\u\\ 2 + C. (3.27) 

When D a = dt, the estimates about I4-I9 can be shown as follows. By using (|3.17p . (|3.21[) 
and the Cauchy inequality, we have 

|J 4 | < C\\Xp t \\ (HU|AVp|| + ||A^||oo||V • n||) < Ck\ (\\XVp\\ 2 + ||Ap t || 2 + ||V • u|| 2 ) , 

(3.28) 

\h\ < C(||Ant||||Vn||oo + HAnlloollVntlDlliitll 

< CdlAntllUVnlla + IIAnllallVntlDKII 

< C||w t || 2 + C||V 2 ri||l + C, (3.29) 



\h\ < C\\v t • VulHkll < C|HU|Vu|||H| < Ciff (||Vu|| 2 + |K|| 2 ), 



(3.30) 



I/7 1 < C"A 2 



Kll < Cll^tlUHAVHIKII < CKi (\\XVp\\ 2 + |h|| 2 ) ,(3.31) 



\h\ < C 



Au 



+ 



V(V-u) 



< OX' 1 (||A6||oo||A«|| + ||A6||oo||V(V • u 



< aA-^KiHi^ + n^ir^, 



(3.32) 



|/ 9 | < c 



- (An ■ Vn) 

Ut 



< CA^IlAetlloollAnlHIVnlUII^II < CA- 2 ^i|K|| 2 + C. 

(3.33) 

Thirdly, by using Lemma 13.14 we give the estimates about 14-/9 when D a = ViV 13 for 
< s- 1, s > 3. 

\h\ < c [Hv^Hooiivvil + livpiuiiv^n + ||ve||ool|v— X (V • «)|| + ||v • ullooliv^n] 
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xA||AV(p-l)|| 



(3.34) 



■ [ V 3 (An • Vn) • AV^u 

< C||V 2 V /3 n||||V /3 (An- Vn) || 

< C||V 2 Vv|| (HAnlloollV'nll + ||Vn||oo|| V^Ar^ 



< ^\\Vuf s + C\\V\\\U, 



(3.35) 



where we have used integration by parts and (|3.17p in f)3.35j) . Similarly as (|3.34[) . one 
obtains 



|I 6 | < C||«|| a (||Vt;||oo||V«|| i -i + ||V«||oo||u||.) <CK?\\u 



(3.36) 



On the other hand, noting the fact that 

'p'(0 



v 6 



< C(m\\s-i + IIV^IU) = C (A-^lAVelU-i + A-'HAV^Ii;.!) , (3.37) 



which can be proved by induction, we have 

£p"(6-no 



|/ 7 | < cx 



< 



V a - L Vp\\ + ||Vp|| 



p'(0 



C [||AV£|| 2 ||AV>|| + l|AVp|| 2 (||AVeiU-i + A 1 — ||AVeil2_i)] 



< c Msr* + \^ki){\\\{ p - i)||2 + |H|2), 



(3.38) 



\h\ < C\\u\ 



V _1 Au|| + ||Au|| 



+ 



vl i 

-1 1^2 



llv-Mv-u)!! + liv(v-u) 



vi i 



< C A- 1 ^ + \- s K? ||«|U(||Vu||. + ||V • u\\ 8 ) 



< C (\~ 2 Ki + X-^K!) \\u\\l + |||Vti||2 + ^||V • u\\l, 



(3.39) 



|/ 9 | < C\\u 
< C\\u 



V 



V _1 (An • Vn)|| + ||An • Vn|| 



V s 



|| Veiloodl An||oo II V—^-VnH + HVnlUIV^AnH) + ||An|| 2 ||Vn|| 2 
< C||u||. [A^HAVeihllAnllallVnlU-i + A" 1 1| AV^|| 2 1| Vn|| 2 || AnIU-! 
+ ||An|| 2 ||Vn|| 2 (A-^IAV^IU-i + || AVeH^O ] 



V 
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< C (\- 2 K x + \- 2s K{) \\u\\l + C||V 2 n||2_i s (3.40) 

where we have used (|3.17p in (|3.4Qj) . 

We now tie the estimates ([3TM - (|3~27l) and (^34]l - (l3T40]l together and get 



d x - 

|a|<s 



W -A 2 |V a (p-l)| 2 + £|V a u| 2 



+ V>/ |VV Q n| 2 + V(/ t + ^) / |V a (V-u)| 2 
< C (K{ + K 2 + l) (\\X(p - 1)|| 2 + ||n|| 2 ) + C||V 2 n|| 2 _ 1 . (3.41) 
Recalling the constraints on the initial data (|'2.8p . we have 

||Apo|| 2 + |K +U || 2 < {2\~ 2 5l + IKH 2 ) . (3.42) 
Then by (j23j) . (f3~T8j) . (j3ld]) . (pTHZj) and the Gronwall's inequality, we get 
l|A(p-l)||^ + ||n|| 2 < e^*^ (V 2 <^^ 

provided that Tq < T3 = min |t2, j^| + ^ 2+1 1. Furthermore, we integrate (|3.4ip over [0, t] 
and get 

p I ||Vu|| 2 + (« + /i) /"||V (3.44) 
jo jo 

Finally, we estimate Z^-ig when D a = ViS/^dt with | — y | < s — 2. For Ig, by virtue of 

(|3.43p and the Cauchy inequality, we are led to 

\h\ < C\\u t \\ s -i (||V«||oo||V'- 2 V«t|| + livntiiooliv 5 - 1 ^!! + iKlloollv'-Vuii 

+ ||V«|| 00 ||V s -Sll) 
< C|M| s _i fiff|K|| s _i+iff \\Vuth + Kl\\u\\ 



< C(K 1 + K 2 + l)\\u t \\ 2 _ 1 + ^\\Vu t \\ 2 + C, (3.45) 



while for the terms I4, -Z5, I7, Ig and Ig, some extra discussions have to be given since 
the methods to estimate these five terms are different between the cases s = 3 and s > 4. 
Indeed, when s > 4, we use Lemma EH the Sobolev embedding # 2 (T W ) ^ L°°(T Ar ), the 
Cauchy inequality and (|3.43p to estimate these terms and get 

\h\ < \\Xpt\\s-i [llVulUllAV^Vptll + IIAVptlUllV^II + II^IUIIAV^Vpll 
+ ||AVp|| 00 ||V s -\|| + ||AVe||oo||V s - 2 (V • ut)\\ + ||V • ntlUIIAV^Cll 
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+ ||A&||oo||V s - 1 (V • u )\\ + ||V • ulUHAV'-'e 



< CHAptll,-! \Kf\\Xpt\\ s -i + Ki\\X(p-l)\\ s + K^\\ Ut \\ s ^+Ki\\u\\ s 

< C{K X + K 2 + 1) {\\\ptt-i + WKp - 1)112 + KIlLi) + c, 



(3.46) 



\h\ 



-v [ V 7 (An t • Vn) • AV 7 u t - v f V 7 (An • Vn<) • AV 7 u t 

< C (IIAntlUHV^Vnll + || Vn|U|| V s - 2 An t ||) ||AV 7 ^|| 
+C(||An|| 00 ||V s - 2 Vn t || + || Vn t |U|V s - 2 An||) ||AV 7 u t || 

< CdlVntHallVnlU-a + HVnllallVntlU-iJIlVutlU-i 
+C(||Vn||3||Vnt||,_2 + ||Vnt|| 2 ||Vn|| a _i) ||V^|| s _i 

< |l|Vn 4 || 2 _i + C||Vni|| 2 _i, 



(3.47) 



where we have used integration by parts and (|3.17|) in (|3.47|) . By using a similar induction 
procedure of (|3.3T[) . we have 



and 



V s 



< c 



< c 



<c(||vei| s - 2 + ||ve||^) 



(3.48) 



p"(0-p'(0 



llv-^tlK civile 



V s 



x /P"(o-p'(a 



< c||6|| s -i + c||6l| s -i(||vei| s -2 + ||v^| 



s-l\ 
s-2) 



(3.49) 



Then we use these two facts to obtain the estimates about I7, Is and L 



\I 7 \ < A 2 ||u t || s _ a 



V 



p'id) 



+ 



IIV^Vpll + ||Vp|| 



V s - 2 V Pt \\ + ||Vpt||c 



i(p'M 



i(p'M 



< C||«f||.-i 



K^+A 2 "^ 2 ||A Pt || s _i + X| A-^+A^ 2 +1 ||A(p-l)|| 



< C + Ki + 1) (IIApiH 2 .! + IKH 2 ^) + C, 



(3.50) 



|J 8 | < CWutWs-! 

+ ||Ali|[oo 



V 



V S - 2 A^|| + ||Au t || c 



1 

I 

V s - 1 V(V-n)|| + ||V(V-u) 



V 



s-1 



+ 



1 



+ 



V 



V S - 2 V(V-^)|| + ||V(V-^)|| C 



1 



V 



s-l 



ys-1 



1 
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< CK|| S _! 



X- 1 ^ + \ i ~ s K l 2 ) (||V«i|| a _i + ||V • utlU-O 



+\~ 1 K 2 2 \- l Kf + A 1 " 5 ^ 2 + 1 ) (||V«|| a + ||V • u 



< g||V«t| 



s-1 



+ 



K + A, 



V • utlG-i + C (K{ + A 2 S + 1) \\utt_! + C (||Vu|| 2 + ||V • u\\ 2 s ) . 



(3.51) 



|J 9 | < C||D^|| 



1 



||V' -i (An- Vn)|| + ||An • Vn|| c 



+ 



||V s - i (An-Vn) t || + ||(An-Vn) 



1 1| oo 



V 



V 



s-l 



s-l 



< c\\v? 



u t 



IIAnHoollV^Vnll + 



1 



+||An||oo||Vn| 



V 



s-l 



+ 



+ 



+ 



V 



||V7l||oa||V*- 2 Ant|| + 



V 



HVnlUHV^Anll 
IIAntHoollV— 2 Vn|| 



V 



IIAnlloollV'^VntH 



+ ||An t ||oo||Vn 



llV^IUHV^Anll + IIAnlloollV^II 

r 



v 



s-l 



V 



s-l 



< C\\ut\\ 



s-l 



A" 1 ( Kl + iff ) + ( A^iff + 1 ) ( \~ l K* + \ 1 ~ S K 1 2 



I" 1 



A- 1 ^ 2 + X S K 1 2 1 (||An t || 2 + ||An t || s _ 2 ) + A K 2 2 HAnlU 
< C {K{ + if| + 1) IKHJU + C (|| VntH^! + ||V 2 n||;U + l) . 



(3.52) 



On the other hand, when s = 3, more refined estimates are needed. We here only deal 
with the terms with \<j\ = s — 2 = 1 {D a = ViVjdt) for simplicity, and the case I7I = can 
be estimated similarly (actually more easily). We use Lemma 13.14 the Sobolev embedding 
H l (J N ) <-> L 4 (T N ), H 2 (T N ) ^ L°°(T N ) and the Cauchy inequality to give 

\h\ < A||AV<VjPt|| [||ViVj(u • V Pt ) - v ■ WiVjptW + HViVj^ • Vp)|| 
+ ||V i V j (£V • u t ) - £V • ViVjUtW + • v«)||] 

< A||AV 2 pt|l [HViVj-u • Vpt + ■ VV iP t + V iV ■ VV jPt \\ + ||V 2 (i> t • Vp)|| 
+||ViV^ • (V • u t ) + V£ • Vi(V • u t ) + V,e • Vj(V • «t)|| + ||V 2 (& • Vu)||] 

< CA||AV 2 p<|| [||V 2 ?;|| L 4||Vp t || L 4 + HVHUIVVH + \\v t ||oo||V 3 p|| + ||Vp|U|V 2 ^|| 
+ ||V 2 e|| L 4||V • u t \\ L 4 + ||V£||oo||V(V • u t )\\ + ||6||oo||V 3 n|| + ||Vu|U|V 2 6||] 

< CA||AV 2 Pt || [||V\||i||V ft ||i + ||Vt;|| 2 ||V 2 p t || + ||^|| 2 ||V 3 p|| + || Vp|| 2 ||V 2 ^|| 
+ ||V 2 e||i||V • ut||i + 2||V£|| 2 ||V(V • u t )\\ + ||6l| 2 ||V 3 n|| + ||V U || 2 ||V 2 &||] 
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< Ciff||Ap t ||2 + cif|||A Pt || 2 + c^f||Av 2 pt||IKIl2 + c'iir|||AV 1 

< c(k x + k 2 + 1) (|| Apt ill + Kill) + a 



u \\3 



(3.53) 



Noting that 



• / [ViVj (An* • Vn) + V«Vj- (An • Vn t )] • ViVj-u* 
-v j [Vj (An t • Vn) + Vj (An • Vn t )\ ■ AVjU t 

-v j (VjAnt ■ Vn + An t • VVj-n + VjAn • Vn t + An • VVj-nj) • AVjU t , (3.54) 



then we obtain 



|/ 5 | < C (|| VAn t ||||Vn|| 00 + 1^^111^^1100 + || AVn|| || VntlU + IIAnlUUVVll) ||AV«t|| 
< C (||VAn t ||||Vn|| 2 + ||An t || ||V 2 n|| 2 + ||AVn|| ||Vn t || 2 + || An|| 2 || V 2 n t ||) ||AVu t || 

(3.55) 



< ^l|V^||i + C(||Vn t || 2 + ||V 2 n|| 2 ). 



Clearly, 
h = 



y / e 

I7.1 + 17,2, 



P'(£) \ P'(£) 
v * v j ( — ~ — ^VViV.-pt 



Vp 



VjVjUt 
(3.56) 



then we get 



M < ca 2 ||vV|| (live HLllvptii + l|v 2 ^|| L 4||Vp f || L 4 + iiveiioollvVii) 

< C||V 2 « t || (||AVe||i||AVpt|| + HAV^IIillAVptllx + ||AV^|| 2 ||AV 2 pt||) 

< c(k 1 + k;A (||Apt|| 2 + |K|| 2 ), 



(3.57) 



where we have used the fact that 



V,V, I ^-^Vpt) - ^-VViVjPt 



1 ^ e ' r v e 

2P"(0 2P'(0 



+ 



e e 2 
e 2 



+ 



v,e • v,e • v P4 + g — ^v<v^ • v Pt 



• Wipt + ^ — — — • vv ^- ( 3 - 58 ) 



We infer from jgttSj) and ([336]) that 

'no 



|/ 7 ,2| < CX 2 \\V 2 u t 



IIVVII + l|Vp|| 



V 



2 , P'(0 



< CK||2||A(p - 1)|| 8 [IIA6II2 + HA&H2 (A^HAVeil + A- 2 ||AV^|| 2 )] 



< c (Kf + k% + 1) \\u t \\ 2 2 + a 



(3.59) 



In conclusion, we have 



|/ 7 | < C (Kf + K\ + 1) (\\X Pt f 2 + \\u t f 2 ) + a 



(3.60) 



Next we come to give the estimate of Ig. 
\h\ < C\\V 2 u t \\ 



+ 



V 2 



+ 



+ 



V 2 



V (V • ut) 



V(V • u) 



V i V i V(V • u t ) 



+ 



V 2 



|) V(V-u) 



< C\\V 2 u t \ 



||VC||L||A^|| + ||V 2 C|| i4 ||A^|| i4 + HVeiUlVAwtll 



+ 



1 

I A 



|V 2 An|| + ||An| 



V 2 



1 

I A 



+ l|veil 2 ollv(v-^)|| 



+ l|V 2 e|| L 4||V(V • ut)\\ L A + ||V£||oo||V 2 (V • u t )\\ + 



+||V(V-n) 



1 



l|V 3 (V ■«) 



v 2 (- 



< C\- l \\V 2 u t \\ {X- l m - 1)111 (HVtitll! + ||V ■ utlli) 

+ ||A(£ - 1)|| 3 (HVutUa + ||V • n t || 2 ) + ||A&|| 2 (||V«|| 3 + ||V • u|| 3 ) 
(||Vn|| 3 + ||V • n|| 3 ) [||A6|| 2 + \\Kth (A- l [|AVe|| + A~ 2 1 1 AX7^ 1 1 ] } 



< |l|V^|| 2 + ^p||V • n t || 2 + C (Kf + Ki + 1) |M| 2 + C (||Vn|| 2 + ||V • n|| 2 ) , (3.61) 
where we have used the fact that 



K + fi. 



^3 VifVjfAut - ^ViV^ViAnt - ^ V^V.Aut - ^V^VjAut, (3.62) 



and similar calculations on the term VjVj 



|V(V-«t)l - |V(ViVjV • « t ). 



For J9, noting that 



e ViV,- 



~(Vn • An) t 



-ViVj(An • Vn) t + VjVj 



^9,1 + -^9,: 



- (An • Vn) 

Ut 



ViVjUt 
(3.63) 



one obtains from (|3.49p that 



1*9,1 1 < C 



V 2 



Ut 



(Vn • An) 



|V 2 n 4 | 
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< c 



||V 2 (Vn • An) || + ||Vra||oo||An| 



|V 2 n t | 



< OX' 1 {||A&|| 2 ||Vn|U|V 2 An|| + ||Ae* II II A^IU || V 3 «|| 

+ ||Vn|| 2 ||An|| 2 [||A£ t || 2 + ||A&|| 2 (A^AV^I + A- 2 ||AV£|| 2 )] } ||V 2 u t | 

< C (Kl + Kl + 1) Kill + C||Vn||§, 



(3.64) 



and 



|/9,2 1 < C(||Ve|| 2 ||(A«-Vn) 4 || + ||V 2 e!|L4||(An-Vn) t || i 4 
+ ||V£|U|V(An-Vn) t ||)||V 2 u t || 

< C (A-^IAVeilillAntllUVnlloo + A~ 2 1| AV^||1 1| AtiIU || Vri t || 
+A- 1 ||AV 2 e||i||An t || 1 ||Vn|| 00 + A" 1 )! AV 2 ^|| HAnH^ || Vri^li 
+A- 1 ||AVe|| 2 ||VAn t ||||Vn|| 00 + A -1 ||AV£||2|| An t || ||V 2 n||oo 
+A- 1 ||AVe|| 2 ||VAn||||Vn < || 00 + A^HAVCIHIAnlUIV 2 ^ 

< C^i+^i||Vn||3 + ^||Vn|| 3 + ^||Vn t || 2 ^ ||« t || 2 

< C (Kf + K^ \\u t \\ 2 2 + C (\\Vng + \\Vn t g) . 

In conclusion, we have 

|7 9 | < C {Kl + Kl + 1) ||u t || 2 + C (||Vn|| 2 + ||Vn t || 2 ) . 



u t 2 



(3.65) 



(3.66) 



Now combining the estimates (pT28|H|333] ) . (l3T45l) - (^47j ). (^50]) - (l3T53| ). (|335jh IpOpl . 
()3.6ip and f|3.66j) together, we have 



d \ -1 

(It ^ I n 

\/3\<s-l J1 



^|AV^ t | 2 + e|V^ f | 2 



+ V nf |vv^| 2 + V + |v-v^p 

||S|<--1 ^ |/3|< S -1 ^ 

< c(#f + jrj + i)(||A ft || 2 _ 1 + |K|| 2 _ 1 ) 

+C (||Vu|| 2 + || V • n|| 2 + IIV^H 2 ^ + \\VntWU + l) . 



(3.67) 



Recalling the constrains on the initial data and (|1.38p 9. we have 

||AV^(*,0)|| 2 + ||vMx,0)|| 2 
< C { || A[C«o + ul) • V]p A ||Li + WKpI + 1)V • H A || 2 _i + l|A 2 Vp A || 2 
+ II IK +A) ■ V](«o +^o)ll'-i + ||A(«q + n A ) + V(V • II^II 2 .! 
+ IIA^o + n^-VK + n^ll 2 .!} 
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s-1 



< C. (3.68) 
Then by the Gronwall's inequality and the fact following by ipTIgjl . (ET221) . flSHD that 



f {\\Vuf s + ||V • u\\ 2 s + ||V 2 n||lU + llVntH^,) < C, 
Jo 



(3.69) 



we get 



(3.70) 



(IIMULi + KIlLi) + M / ||V«t||2_i + (« + M) / IIV • ut||2_i < c, 

jo jo 

provided that To is small enough such that Tq < T4 = min j T3 , K s + 1 K s +1 j • 

It remains to show the first inequality of (|3.1|) . It suffices to show ||A(p — l)|| s + \\u — 

xio || s + ||Vn — Vno||s-i < Cq 1 5 by Sobolev's inequality, where cq is the Sobolev constant. 

Let ~p = p — 1, u = u — uq, Vn = Vn — Vuq. We proceed as (|3.4p . ()3.9[) and (|3.10p . and 

then use fl277D, (JZSD and (l3TTTjl to give 



^nll 2 



I 2 ||VV^n|| 2 + Yl ll y2v/ 

|/3|<s-l |/3|<«-1 
< CK^Vnf^ + CK 2 (||n|| 2 + ||n || 2 ) + C {\\{v ■ V)n || 2 _! + ||An || 2 _i + II |Vm| 2 n || 2 _i) 



Is-] 

< CK X || Vn|| 2 _! + C(Kf + l), (3.71) 

where we have dealt with the term such as J Tn VjV^(u ■ Vno) • VjV^An by integration by 
parts and the Cauchy inequality as follows: 

/ V i V' 3 (u • Vn ) • ViV^An = - ( V?(v ■ Vn ) • AV^An 

Jfn Jjn 

< luAV^nf + CHV^^Vno)!! 2 , (3.72) 
and the terms f Jn V i V /3 (An ) ■ V^V^An and J T „ ViV^(|Vm| 2 n ) • V^V^An were treated 



by the same method showen in (|3.72p . Since (|2.8p implies that || Vn(s, 0) \\j < A <5 2 , we 



conclude 



|Vn|| 2 _! < [\- 2 5l + C(K\ + 1)T ] e c ^ lT ° < c^, 



(3.73) 



provided that A 1 and 7b(< T4) are both sufficiently small such that (I3.73P holds. And 
then integrating (13.7ip over [0, t], we have 



/ ||V 2 n|| 2 dr < \- 2 5l + CK^Tq + C(K\ + 1)T . 
Jo 



(3.74) 



Similarly as (|3.4ip and (|3,7ip . we have 

|E f^|Av"p| 2 + eiw 2 



a|<s 
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+ Vu/ |VV a u| 2 + V (« + //) f |V a (V 



| ct <s 

< C{K{ + K 2 + l 



+C 



+c 



^"s + INI') + C||V 2 n||Li + C||(t; • V)u f s + C||V«o| 

2 

2 



V s Au 



V 6 



V (V • u ) 



+ C||V • u | 

2 



V S V(V • u ) 



< C(K[ + K 2 + 1) 



+ \\u\\ 2 s ) + C||V 2 n|| 2 _ 1 +CK X + C. 



(3.75) 



Then following by (|2.8p , (|3.74p , (|3.75p , the Gronwall's inequality and the Cauchy inequality, 
we get 



< e 



C(K[+K 2 +i)T 



2\- 2 5^ + C I ||Vn|| 2 + C(K X + l)r 
Jo 



< Ce c ( K i +K * +1 ) To [\- 2 5 2 + CAWo + (AT 2 + l) T ] 

< c^d, (3.76) 

provided that Tq(< T4) and A -1 are both small enough such that (|3.73p and (|3.76p hold. 
As a conclusion, we have the following lemma. 

Lemma 3.2 Under the assumptions of Theorem \2.1l and soppose that Bj< (Uq) is defined 
by \3. 1\) and A : V —¥ U is defined by the system jl3.fy) . Then there exist constants To, 5, 
K\ and K 2 independent of X such that A maps (Uq) into itself. 

Now we plan to show that A is a contractive map in some sense. In the proofs of this part 
and the following lemmas in this section, C denotes the constant depending on the initial 
data, the domain, N, s and the viscosity coefficients p,, k, v and 9. Firstly, let U = A(V), 
U = A(V), where V, V E Bj, . Then by the definition, we have 

(p - p) t + (v • V)(p - p) + £V •(«-«) + [(v - v) ■ V]p + (£ - f)V ■ u = 0, 
(u - u) t + (v-V)(u-u) + [(v - v) • V]u + A 2 ^V(p - p) + A 2 - ^) Vp 



= f A(u - u) + (J - |J Au + ^pV[V • (u - u)] + - V(V ■ 

-|[A(n - n) ■ Vn + An • V(n - n)] - (| - |J An ■ Vn, 
(n — n)t + (u • V)(n — n) + [{v — v) ■ V]n 

= 8A(n — n) + 0|Vm| 2 (n — n) + 9[(Vra — Vm) : (Vm + Vm)]n. 
Multiplying ()3.77j) ^ by (n — n) and A(n — n) respectively, we have 



(3.77) 



ld_ 

2 jx n 



n — n\ + 



72 — Tl\ 

(v ■ V) h [(f — v) • V]n • (n — n) 



2fi 



-6 [ \V(n-n)\ 2 + 9 f [Vm| 2 |n-n| 2 

Jfn Jfn 

+9 I [(Vm + Vm) : (Vm — Vm)]n • (n — n) 



(3.78) 



and 



--— f |V(n-n)| 2 + [ (y- V)(ra-n) • A(n-n)+ / [(v - v)V]n • A(n - ra) 

2 »' Jfn Jfn Jfn 

= 01 \A(n-n)\ 2 + 6 \Vm\ 2 (n - n) ■ A(n - n) 

+6> / [(Vm - Vm) : (Vm + Vm)]H ■ A(ti - H). (3.79) 

Then multiplying (|3.T7[) i and ()3.77p 9 by A 2 P ^ (p — p) and £(it — u) respectively, we have 
A 2 d f P'(£), ~ |2 A 2 f n 2fP'(0\ 



+A^ 



p'(0 



T" 



c ( P -p)(v- V)(p -p) + \ 2 P'{Z){ P - p)V -(u-u) 



+A 2 / ^l(p-p)[(v-v)-V]p + X 2 



no 



(p-p)(^-OV-S = 0, (3.80) 



and 



it — u 



\~u I £.\ u -u\ 2 -\f Zt\u-u\ 2 +( f(u-V) 



+ / - v) ■ V]u ■ (u-u) + X / P'(£)[(u - u) • V](p - p) 

Jfn Jfn 

-li [ |V(u -u)\ 2 + I e f 7 - -1 AS- (u - u) - ( K + fi) f |V • (u - u)\ 



+ ^e(^-^][(«-8)-V](V.«) 



/ [Vn • A(n - n) + V(n - n) • An] ■ (u-u) + £ ( ^ - ^ ) [(u - u) ■ V]n • An. 

JT™ ./T™ 



V V 



(3.81) 



In conclusion, using H3.78n (j3.8ip and preceding similarly as before, we obtain 



ld_ 

2 (it ./jn 



p'(£) 



X(p — p)\ 2 + £|it — u| 2 + |n — n| 2 + |V(n — n) 



+u / |V(u-u)| 2 + (tt + u) / \V ■ (u-u)\ 2 + 6 [ \V(n-n)\ 2 +6 [ |A(n-n)| 2 

Jfn Jfn Jfn Jfn 

< C (\\X(p - p)\\ 2 + \\u - n|| 2 + ||n - n\\\ + ||A(£ - 6II 2 + ||« - £|| 2 + ||V(m - m)|| 2 ) 

ft 

(3.82) 



+fl|V(n-S)|| 2 + ^||A(n-n)|| 2 , 
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where we have used the following two estimates different from before and the same discus- 
sions also have been used to deal with the last two terms of the R.H.S. of ()3.81f) . 



< C||£ — £|| || A2||£4||it — u\\ L 4 

< C\\t- f||||A2||i(||it - u|| + ||V(u - u) 

< c||e-£]|(||ti-u|| + ||v(u-u)||) 



< 



l l i 



V(u -u)\\ 2 + C(\- L \\\(£ - Oil' + \W ~ «lr),(3.83) 



and 



e(^-^)V(V.S)-(«-2) 



<^iiv(«-s)f +c(\\\(c-or +\w-u\n 



(3.84) 



Then noting that (U — U)(0) = 0, the Gronwall's inequality, and ()3.82f) . we obtain 



sup (||A(p-p)|| 2 + 

0<t<T 

o CT 



2 , II — 1 1 2 , || ~||2\ 

u — u\\ + ||n — n\\i) 



< CT e cl ° sup [\\\(£-0\\ +\\v~v\\ +\\m-m\\( 

0<t<T 



(3.85) 



Moreover, we integrate (|3.82p over [0, t] and get 

\u — u\\1 + 1 1 77. — fi 1 1 







< C (T 2 e cl ° + T ) sup +\\v-v\r + \\rn-rh\\i 

0<t<T 



(3.86) 



If we take T (< T 4 ) small enough such that (|3T73jl . (137761) and C (T$e CT ° + T e CT ° + T ) < 1 
hold, then we can prove the contraction. 

In conclusion, we have the following lemma. 

Lemma 3.3 Under the assumptions of Theorem \2.1\ the maps A : V — > U is a contraction 
in the sense that 

r-t 



sup (||A(p - p)|| 2 + ||u - u|| 2 + ||n - n|| 2 ) + / 
<t<T Jo 



0<t<T 

< t sup ( ||A(£ — £) || 2 + \\v — v\\ 2 + \\m — m\\f 

0<t<T ^ 

for some < r < 1, provided that Tq is small enough. 

In order to prove Theorem l2.lt we also need the following lemma. 



\u — u\\1 + \\n — nil 2 ) 



(3.87) 



Lemma 3.4 Consider the incompressible system of liquid crystal |1.^Q[ ) with the initial 
condition j2.7\) for s > 3. Then there exists at most one strong solution. 



Proof. Assume that (ui,m) and (1*2,^2) are two strong solutions of (jl.40p with the same 
initial data (|2.7p . Then we have 

V • (ui — u 2 ) = 0, 

d t (ui - u 2 ) + («i • V)(ni - u 2 ) + [(ui - u 2 ) ■ V]u 2 + V(pi - p 2 ) 
= /iA(m - u 2 ) - vV ■ [V(m - n 2 ) VnJ - i/V • [Vn 2 V(n x - n 2 )] , ( 3 - 88 ) 
<9j(rei - n 2 ) + (m ■ V)(ni - n 2 ) + [(ui - u 2 ) ■ V]n 2 

= 0A(m - n 2 ) + ^iVnil 2 ^! - n 2 ) + 0[(Vni - Vn 2 ) : (Vnj + Vn 2 )]n 2 . 

Multiplying (I3.88P 9 by {u\ — u 2 ) and multiplying (j3.88jU by {n\ — n 2 ) and A(ni — n 2 ) 
respectively, we have 

V(«i - u 2 )|| 2 + 0||V(m - n 2 )|| 2 + 0||A(m - n 2 )|| 2 
< C (||m - u 2 \\ 2 + ||m - n 2 \\\) + |||V(m - u 2 )|| 2 + -||A(m - n 2 )|| 2 . (3.89) 

Noting that ||tti(x, 0) — « 2 (x,0)|| 2 + ||ni(x,0) — n 2 (x,0)|| 2 = 0, and by the Gronwall's 
inequality, we have 

||«l - u 2 \\ 2 + ||m - n 2 || 2 = 0. (3.90) 

This completes the proof of Lemma 13.41 □ 

Before proving Theorem 12.11 we give the following lemma. 

Lemma 3.5 (128^ ) Assume that ICficY are Banach spaces and X r > r ) E. Then the 
following embedding are compact: 

(i) ^ip:<peL q {0,T;X),^eL 1 (0,T;Y)\^^L' 1 {0,T;E), if 1 < q < 00; 

(m) Lp:<peL°°(p,T;X),^ eL r (0,T ; y)} C([0,T\- E), if 1< r < 00. 

Now we are in a position to give the proof of Theorem 12.11 

Proof of Theorem 2.1. For any fixed A, the standard procedure produces a sequence 
{(pijUijTii)}™! satisfying 

dtrf+i + K • V) p} +1 + p A V • u} +1 = 0, 

dtul, + K • V) u A +1 + A 2 ^Vp A +1 = 4A^ +1 + ^1 VV • nf +1 - f (An A +1 . Vn A +1 ) , (3.91) 
dtnf +1 + (u A • V) n A +1 = 6 (An A +1 + |Vn A | 2 n A +1 ) , 
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as well as the uniform estimates 

ft r- 







E s (U x (t)) + / ^l|V«^ + (« + //)||V-^||2 + 6»||V 



n 



A 1 1 2 



<c, 



E a -i{dtU^(t))+ / ^llv^il^i + C/c + ^llv-^liLi + eilv^n^ 

Jo L 



(3.92) 



< C. 



Let pf +1 = p x +1 — p x , u x +1 = u x +1 - u x , and n x +l = n x +1 — n x . In view of Lemma [3731 we 



have 

OO 



OO CO 

All ^ „„ / II^A 



|Pj || < oo, > ! | u 

i=2 i=2 



, E(ll^ A H + / V^fJ <°°> E(|I^ A Hi + / °II A ^II 2 ) <oo.(3.93) 



oo oo 

Let p x = p\ + Yl Pi ') nA = u i + ]C an< ^ nA = n i + S "- A ' then we have 

i=2 i=2 i=2 

p X ^p X , in L°°([0,T ];L 2 ), 

u x ^u x , in I 00 ([0,ro];L 2 )nL 2 ([0,ro];fl- 1 ), 
n A ^n\ in L°°([0, T ]; tf 1 ) n L 2 ([0, T ]; tf 2 ). 
It follows obviously that 

p\ u\n x G L°°([0, T ]; tf s ) PI Lip([0, T ]; If- 1 ) 



satisfy the estimates (|3.92p according to the lower semi-continuity. For any s' € [0, s), by 
the Sobolev interpolation inequalities, we have 

\\(p X , U X ,n X )-(p\ U X ,n X )\\ s , 
< C\\(p x ,u x ,n x ) - (p x ,u x ,n x )\\ e (\\(p x ,u x ,n x )\\ s + ||(/, u\ n^) 1 ' 6 (3.94) 



as i — 7- oo, where we have used lemma [3731 to get 



\\(p x ,u x ,n x ) - (p\u\n x )\\ < \\{p^u x ,n x ) - (^_ 1 ,«}_ 1 ,nJ_ 1 )|| < — . (3.95) 



Hence, (p x ,u x ,n x ) E C([0,T ]; H s '). Besides, with the aid of Lemma [375] and (13.92j) . one 
deduces easily that (p x ,u x ,n x ) is a strong solution of compressible Liquid crystal model 
(|1.38p . Finally, multiplying (|l,38p q by n x , we get an equation for (|n A | 2 — l) as 



i^Y.jj + I (V. V J ^| n A|2 _ ^ = ^A(|n A | 2 -l)+e|Vn A | 2 (|n A | 2 -l) (™G) 



Multiplying (137961) by (|n A | 2 - l), integrating over T n , and then using the Gronwall's in- 



equality and the assumption that |n A (x,0)| 2 = 1, we have 

2 

for all t € [0, T, 



1 


|n A | 2 -l 







0J 



(3.97) 
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which implies \n x \ = 1 a.e. in Qt , then from the regularity of n x , we conclude that 



\n | = 1, in Q To . 



(3.98) 



The uniqueness can be proved by a similar argument as Lemma 13.31 This completes the 
proof of the uniform stability part of Theorem 12. 1L 

Next, we are going to show that (p x ,u x ,n x ) converges to the unique strong solution to 
the corresponding incompressible system (jl.40p as A — >■ oo. To see this, note first that (|2,9p 
implies that p x — > 1 in L°°([0, To]; H s ) n Lip([0, To]; -fP -1 ), and there exists a subsequence 
{(u Xj , n A J )}j of {(u x ,n x )}\ with a limit u and n such that 



(u x i,n x i) (u,n) weakly* in L°°([0, T ]; H s ) n Lip([0, T ]; il 



s-l> 



(^,n A 0^(n,n) in C([0, T ]; # s ') 



(3.99) 



for any < s' < s, where we have used the fact that the embedding H s <^-> H s is compact 

and Lemma 13.51 Now we are to let j — > oo (that is Aj — > oo) in (|1.38p . First of all, 

multiplying (jl.38P i and (|1.38p ^ by two smooth test functions ipi(x,t) and tp3(x,t) with 

compact in t € [0, To] respectively, we have 

rib j- rT /■ 

V • u Aj ipidxdt = / / 

To 



p t ] + (u A i -V) + (p^ - 1 V • -u^ Vi^^, (3.100) 



A, 



+ ( u X J ■ V ) 



At I 2 _„ A r 



^dxdt = 0. 



(3.101) 



/0 JT" 

Then (u, n) satisfies (jl.40P i and (jl.40jh . Let ip(x,t) be a smooth test function of (|1.38p 9 
with compact supports in t £ [0, To] and the divergence free condition V • tfj = 0, we have 

n 

JO JT" 



u^ 3 + ( u A > ■ V ) u A ^ - -^-Au^' 



p X i 



p X i 



V ( V • u Aj 



lb 



H — V • Vn Aj Vn A ^ r-V 

p A i V / p A i 

A,- 



A, 1 2 



■ipdxdt 



rp 3 p'(c) 
-AfVV / — f±d£dxdt = Q, 

*/T n Jl ? 



(3.102) 



then let j — >• oo and get 



u t + (u • V)u - /iA-u + i/V • (Vn Vn) - v V 



|Vn| 



(3.103) 



where P is the L 2 -projection on the divergence zero vector fields. If 



u t + (u- V)u - pAu + zA7 • (Vn Vn) - z^V 



|Vn| 



-Vp 



(3.104) 
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for some p £ L°°([0,T ]; H 3 - 1 ) (1 L 2 ([0,T }; H s ), then we have 
A 2 

-^-V (P(p Aj )) -> Vp, u>ea£%* in L°°([0, T ]; tf s ~ 2 ) n L 2 ([0, T ]; H 8 ' 1 ). (3.105) 

Taking p = p — z ^ V ^ , then (|1.40p 9 follows from (|3.105j) directly. Actually, Lemma 13.41 
ensures that the sequence (p x ,u x , n x ) itself converges as well. Moreover, similarly as (|3.96p - 
(|3.98p . we can prove |n| = 1 in Qt - This completes the proof of Theorem 12.11 □ 

Remark 3.1 It follows from H3.105}) that 

||A 2 Vp A || s _ 2 + f ||A 2 Vp A || 2 _ 1 < C, t e [0,T ]. (3.106) 
Jo 

Remark 3.2 It follows from i fOgj) and IO|) that 

||AV-u A || s _i < ||Ap A || s _ 1 + ||A( U A -V)p A || s _i + ||A( / o A -l)V-u A || s _ 1 <, t e [0,T ]. (3.107) 
4 Global existence of strong solutions 

In this section, we devote ourselves to getting a priori energy estimates (I2.13|) and (|2.14p 
for small initial displacements and small initial data, and then together with Theorem 12. 1[ 
Theorem 12.21 can be proved by a standard procedure. First of all, we assume that |n| = 1 
for all (x,t) € Qt, and 

F s (U(t)) + f [p\\Vu\\ 2 s + (k + /i)||V • + ^llV^ll^i] < 4 (e + A" 2 5 2 ) (4.1) 
Jo 

for i € [0, T], then what we need to do is to prove the following desired estimates: 

F s (U{t)) + f [/i||Vu|| 2 + (k + /z)||V • n|| 2 + 0||V 2 n|| 2 _ 1 ] < 2 (e + \- 2 5 2 ) . (4.2) 
■/ o 

for t € [0, T]. Then (|2. 13j) follows by the standard continuity argument and the fact 
F s (U(0)) <A(e + X~ 2 5 2 ). 

Now we plan to give (|2.14p under the assumptions that \n\ = 1 and ()4.ip . We go back to 
(I3T7D . (^81) . (f3TTTT) f|3TT5]) . (l3T24l) - (^26|) . (l3T28l) - (^33|) and (I3^5l) - (l336|) . replace (£,v,m) by 
(p,u,n), and then we make a similar (just a little different) argument to give the following 
fact: when D a = dt, we obtain 

2 9 

\Nj\ + £ < C (||V 2 n|| 2 + 1) (HAptll 2 + |K|| 2 + K|| 2 ) , (4.3) 
i=i j=i 

32 



where we have used (jl.38P i to give the following fact and then used it to estimate I\: 

, < ||V • (HIU < C(||n|| 2 ||Vp|| 2 + HpIUIIV • «|| 2 ) < CA- 1 + C < C. (4.4) 



On the other hand, when D a = VjV 7 9f with (7! < s — 2, we have 
2 9 



< fnvu t |i2_ 1 + ^iiv.« t |i;_ 1 + 



- + C (e + \- 2 ) 



l|Vn t ||ti 

+C (||Vu|| 2 + ||V • + IIV^H^ + 1) (\\XptWU + \\utWU + IKHti) , (4-5) 



where we have used the fact indicated by (|4.ip that ||Vn||g_ 1 < C(eo + A 2 ) 
Combining (|4.3p and (|4.5p . we have 



I £ i(^« + HvW + |vW) 

|/3|<s-l 1 V r 7 

+ Y (if |VV /3 n 4 | 2 + V {k + h)I |V-V /3 n t | 2 + V (9 / |VV / 



2 



l/3|<s-l R \0\<s-l |/3|<s-l 

< C (||Vn|| 2 + ||V • n|| 2 + IIV^H 2 ^ + 1) {\\\pt\\U + \\u t \\U + \\ntfs-l) , (4-6) 

provided that Sq and A -1 are both small enough. 

Then by the GronwalPs inequality, (f2~5]) . ([3720]), ([BTggj) and (}4~f|) . we get 

IIA^H 2 .! + ||n t || 2 _i + IKIIi-i + f [/^HVn t || 2 _i + (« + M)I|V ■ «*H2— 1 + ^V^H^i] 

JO 

< C(l + i)expCt < CexpCt, t £ [0,T]. (4.7) 

And thus we can conclude that Remark [3 . 1 1 and [3721 hold for t G [0, T] by the same discussion 
shown in Section 3. 

Secondly, we shall turn to prove (|4.2j) . We go back to (13.41) . (|3.9p and (I3.10p . replace 
(v,m) by (u,n), and then make a different argument by noting that \n\ = 1 to give 

= f V p (u- Vn) • / (|Vn| 2 n) • AV' 3 



< C||V 2 V /3 n|| ( ||V^(u • Vn)|| + ||V^ (|Vn| 2 



n 



2 n) II 2 



< |||V 2 V /3 rz|| 2 + CHV^u • Vn)|| 2 + C\\V P (|Vn| 

< |||V 2 V^|| 2 + C (||n||Ll|V^|| 2 + HVnllLH V-^f) 
+C (llVnll^HV^nll 2 + llnllLllVnllLUV^II 2 ) 
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< ^iv^nll 2 + C (|M| 2 ||V s n|| 2 + || Vn||l|| V'-^H 2 ) 
+C (||Vn||f llV^nf + ||Vn||l|| V'n|| 2 ) 

< |||V 2 V^|| 2 + C{e + A- 2 J 2 ) ||V s+1 n|| 2 , 

where we have used Poincare inequality in the last step since 



(4.8) 



V k n = 0, k > 1. 



(4.9) 



We conclude from (14.81) that 



(4.10) 

|/9|<s-l |£|<s-l 

provided that £o an d A -1 are both small enough. Then by integrating (I4.10P over [0, t], we 
have 

HVnH 2 ^ + 6 f || V 2 n|| 2 _i < HVnoll 2 ., + A~ 2 5 2 < e + A-% 2 . (4.11) 
J o 

Thirdly, recalling (|3.23p and replacing (£,t>) by (p,u), it yields 



1 d 

2 (it ./x n 



A2^|D«(p-l)[2 + p[DV 



+fi [ \VD a u\ 2 + (k + n) f \V-D a u a 
= Ji + J 2 + h + Ji + -h + J% + Ji + Js, 

where there are some slight changes from Ij shown as follows: 
A 2 f 2P'{p)-P"(p)p^ 



(4.12) 



Ji = — 

J 2 = A 



2 jt 

2 



P 



-V -u\D a {p-l)Y 



■h 
J* 

■h 



-A 



/ P"(p)D a {p-l)D a u-Vp 
P'(p) 



T n /J 



L> Q (p - l){[D a (u • Vp) - « • VDV] + [£> a (pV • u) - pV • D a !i]}, 



-i/ / D a (An ■ Vn) • D a u, 

- [ p [D a (u ■ Vu) - u ■ VD a u] ■ D a u, 



A" / p 



p 



+ / PS D 

It 



D a ( - —AD C 

P J P 

K + p, 



U 



P 



-V(V-u) 



« + A* 



VL> Q (V ■ u) \ ■ D a u, 
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P \D C 



P 



-(An • Vn) 



P 



-D a (An • Vn) } ■ D a u, 



since we have used (|1.38P i to give 

\D% P -l)\ 2 d t (^-)+^ [ |^(p-l)| 2 V 

Z Jjn VP/ Z il" 

\D a {p-l)\ 2 

\ P /J 
fP'(p)\' _ 

Pt + 



P'(P). 



pt+v-[ n 



P J VP 
P'{P)\ , fP'(p)^' 



2 J ipn 

-/ 

2 Jt« L P VP 

^ r 2P'(p)-p"( P )p 

2 j T n p 



+ 



P J VP 
^'(P) (P'{P)\ 



Vp • n + I p\7 ■ u 



P VP 
P'(p) (P'{P)\ 



P 



P 

p V -u\D a (p-l)\ 2 
V-u\D a (p-l)\ 2 



V-u\D a (p- 1)| 2 . 



(4.13) 



Now we give the dispersive estimates about </& for k = 1, 2, . . . , 8, when D a = V a with 
\a\ < s, s > 3. 



|Ji| < CA-iAV-nlloollAV"^-!)!! 2 < CA^HAV^p - 1)|| 2 , 



(4.14) 



where we have used Remark 13.21 



\J 2 \ < CA||Vp|| L 4||V Q n|| L 4||AV Q 09-l)|| 

< CAllVplK (||V a n|| + ||VV a n||) ||AV a (p - 1)|| 

< CA- 1 ||A 2 V /9 || S _ 2 (||V Q n|| + ||VV a n||) ||AV a (p - 1)|| 

< ^||VV Q n|| 2 + CXA- 1 + A- 2 )(||V a n|| 2 + ||AV a (p - 1)|| 2 ) 



(4.15) 



where we have used the Sobolev embedding H 1 (T N ) L 4 (T N ) and Remark 13.21 Then we 
use Lemma [3TT~I to give 



|J 3 | < CA||AV a (p- 1)|| (||V a (n- Vp) -u- VV a p\\ + ||V a (pV-n) - pV ■ V a n||) 

< CA||AV a (p - l)||(||V«|U||Vp|| s _i + ||Vp|U||V s -V|| + HVplloollV^V • n|| 
+ ||V-n|| 00 ||V s p||) 

< CA- 1 ||AV a (p-l)||||Vn|| s _ 1 ||A 2 Vp|| s _ 1 

< CA- 1 ||AV a ( /3 -l)|| 2 + C7A- 1 ||A 2 Vp|| 2 _ 1 , (4.16) 



|J 5 | < C||V a u||||V Q (n- Vn) -n- VV a n|| < CllV^ljllVnlloollV'ul 
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< C £ 2 +A-% l|Vu||ti 



|J 6 | < C\ 2 \\V a u\\ 
< C\ 2 \\V a u\\ 



v 



p 
p 



p'(p) 



VV°> 



|V s p|| + l|Vp||c 



V 



< C\ 2 \\V a u\\\\Vp\\ 2 \\V s p\\ 

< CA- 1 ||V Q n||||A 2 V /0 || s _i||AV s p|| 

< CA- 1 ||V Q n|| 2 + CA- 1 j|A 2 Vp||2_ 1 , 



(4.17) 



(4.18) 



|J 7 | < C\\V a u\ 
< C\\V a u\ 



V a I ^-Au 
P 



V 



^AV a u 
P 

7S-1 



V 



HV^Aull + ||A«| 



a l K + P 

V P 



V(V-n) 



V(V • u) 



+ 



V 



K + P 
P 



l|V(V-«)ll + 



V(V-u)||oo||V 



p 

s l'^jfp 
p 



< CA-^lV^IHlAVplU-illV^II,-! 

< cx- 1 (|| V«|| 2 + ||Vi|| 2 ) , 



(4.19) 



|J 8 | < C||V Q u| 
< C||V a u| 



p 



■(An • Vn) 



P 



-V a (An- Vn) 



VI* 

P 



V _1 (An • Vn)|| + HAn^H Vn|| 



VI ^ 



< C\\V a u\\ [HVHbdlAnlUllV^II + ||Vn|UI V^AnH) + ||An||.||Vn|| 8 ||Vp|| a _i] 

< CX' 1 (\\u\\ 2 s + || V'nWU) . (4.20) 

Particularly, we have to make some argument about the estimate of J4. When |a| = 0, the 
similar estimate shown in (|3.27p can be obtained as 



\J4\<C\\ U \\\\Vn\U\An\\ 2 <C\\ U \\ 2 + C(eo + X- 2 5 2 ) ||V 2 n|| 2 . 



(4.21) 



On the other hand, when V a = VjV 3 , |/3| < s — 1, we give the following estimate similar 

as fllTMD 



[ V /3 (An • Vn) • AV^ 



< ^||V 2 Vu|| 2 + C|| V(An • Vn)|| 5 



< ^||V 2 V /3 n|| 2 + C||An|| 2 ||Vn|| 2 + C||Vn|| 2 ,||V s+1 



ii\ 



< ^||V 2 V^|| 2 + C(e + A- 2 5 2 )||V 



2^112 



n 



Is— 1' 



(4.22) 



In conclusion, we combine (|4.14p - (|4.22p to give 



1 d x - 

2dt ,V 

\a\<s 



P 



s-1 



/ |vv%| 2 + ( K + / ,) V / |v-v c 

\et\<s |ck| <s 

< ~ l|Vn|| 2 + C (e\ + A" 1 ) \\Vuf s + C ( eo + A" 2 ) ||V 2 n|| 2 _ 

CA- 1 ||A 2 Vp|| 2 _i + OX' 1 {\\u\\ 2 + ||A(p - 1)|| 2 ) . (4.23) 

then we can use the Gronwall's inequality, Remark 13.11 and ([4. lip to give (|4.2p . provided 
that £q and A -1 are both small enough. Finally, \n\ = 1 in Qt can be proved by repeating 
the procedure shown in ()3.96j) ()3.98[) . 

5 The convergence rates about u x and n A when A — > 00 

In this section, we will prove Theorem 12.31 by the modulated energy method with the help 
of uniform estimates (|2.9p . 

Proof of Theorem 2.3. First of all, let's rewrite fll .381) as follows: 
pi + div(p A -u A ) = 0, 

(p x u x ) t + V • (p x u x <g> u x ) + X 2 VP(p x ) = pAu x + (k + p)\/divu x - v(An x ■ Vn A ), (5-1) 
n$ + (u x ■ V)n x = 6>(An A + \Vn x \ 2 n x ). 

Multiplying (|5.ip 9 by u x and then integrating over T n , we have 

1/ (Ip X W X \ 2 + \ 2 Q(P X )) + P [ \Vu x \ 2 + ( K + p) [ |V V| 2 

= -v (u x ■ V)n A • An A , (5.2) 

where Q(p) = p ff Multiplying (|5.ip ^ by v(An x + |Vn A | 2 n A ), and then integrating 

over and noting that |n A | = 1, we have 

^- [ -\Vn x \ 2 + ve[ |An A + |Vn A |V| 2 = v [ (u X ■ V)n A • An A . (5.3) 

dt Jjn 2 Jjn Jfn 

Combining (|5.2p and (|5.3p . we obtain 

al Jfn \ z 1 J JY n 

+ {K + P)f \V-u x \ 2 + ve( |An A + |Vn A |V| 2 = 0. (5.4) 
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Let 



n\x,t) = X 2 Q(p x )-P(l)(p x -1) 



(5.5) 



then we have / II A (x,0) < CA 2 from the Taylor series and (|2.8[) . Integrating (I5.4p over 
[0, t] and using the conservation of mass, we have the basic energy estimates as follows: 

/ (~pV| 2 + ~|VnY + n\z,t)W f I |Vu A | 2 + ( K + / u) f I |V-u A | 2 

Jf n \2 2 J Jo JT n JO JT™ 

+v9 f f \An x + \Vn x \ 2 n x \ 2 = ( f^p A |u A | 2 + ^|Vn A | 2 + U x (x, 0)) . (5.6) 
Jo Jf n Ji n \ z 1 J 

Similarly, from (jl.38p . we get the basic energy law for the incompressible system as 

I (-\u\ 2 + -|Vn| 2 ) + p f I \Vu\ 2 + vo! f |An+|Vn[ 2 n| 2 
Jj n \2 2 J Jo Jj n Jo il" 

(5.7) 

Secondly, multiplying ()5.ip 9 by u and then integrating over T n , we have 

p x u x [(u • V)u + Vp - pAu + uV ■ (Vn Vn)] 



iK| 2 + ^|Vn | 2 ). 

qpn \ Z Z 



p x u x u + 



JT n 




(p X u x ®u X ) -\7u + p [ [ Vu A •Vtt + i/ / / (n • V)n A • An A 
'0 -/T 11 JO Jf n Jo Jf n 

PoU X u . 



(5.8) 



Then multiplying (|5.ip ^ by vn, we obtain 




/ n x ■ n — v I n ■ uq + v l l n 

It™ Ji n Jo Jt 

+v I I (u x -V)n A -n = v9 [ [ ^- x 
Jo Jt™ Jo Jj n 



X [(u ■ V)n - 9 (An + |Vn[ 2 n)] 
(An A + |Vn A | 2 n A ) • n. (5.9) 



Similarly, multiplying ()5.ip ^ by vAn, we have 

-v [ Vn x -Vn + vf [ Vn A • {-V [(u ■ V)n] + 9 [VAn + V (|Vn| 2 n)] } 



+v [ [ (u x ■ V)n x ■ An - v6 [ [ (An x + \Vn A \ z n A ) ■ An 
Jo ii" Jo Jf n 



,A|2„A>> 



-1/ 



Vn A • Vn . 



(5.10) 



In conclusion, combining (|5.6p - (|5.10p . and noting the fact that \n x \ = \n\ = 1, we have 

/ (\\^u X -u\ 2 + Un x -n\ 2 + UVn x -Vn\ 2 + Ii x (x,t)\ + p ( I |Vu A - Vu\< 
Jjn \2 2 2 J Jo JT n 



+ (K + p) I f |V-n A | 2 + ^ 

Jo Jt™ Jo JJ n 



(An A + \Vn x \ 2 n x ) - (An + \S7n\ z n) 



lis 



' Po u o - u o 



+ - ^o| 2 + !^|Vn A - Vn | 2 + n A (x,0) 



-H I I Vu x -Vu+ [ \/^(\/^ - l)u A • u - I J,^ ! J/>fi - I I t,<) ■ » 

JO JT" JT" JT™ V / 

+ [ [ p X u x ■ [(u • V)u + Vp - nAn + vV • (Vn Vn)] 
Jo JT™ 

n(p X u x (& u X ) ■ Vu + v [ [ (u x ■ V)n A - n + u [ [ (u ■ V)n • n x 
Jo JT" JO JT" 

-vO [ [ {An + \Vn\ 2 n)-n x -v6 ( [ (An x + |Vn A | 2 n A ) ■ n 
Jo jt™ Jo jx™ 

,A A „A ,, / / / A V7\„A 



,A 



+i> / / (u • V)n A • An A - zW / (u • V)n A ■ An-u / (u • V)n • An" 

JT" io JT n Jo il" 

vO f f (An+\Vn\ 2 n) -\Vn x \ 2 n x - v9 [ f (An x + \Vn x \ 2 n x ) ■ \Vn\ 2 n 
Jo jt™ Jo ip 

+ ^|n A -n | 2 + ^|Vn A - Vn | 2 +n A (x,0) J + £ i? A (t),(5.11) 

/ 1<«<8 



5 



where 



#i (*) = / \/p A (\/p A - 1)« A • u - I Jp x (J p x - 1 j • n , 

JT™ JT n V / 

R x {t)= f [ p X u X • [(« • V)u] - / / (pV u A ) • Vn, 
70 JT" JO JT™ 

JO JT n 

/"* /" /"* f 

A/+\ — ,, / / V7„,A V7„, ,, / / „A„,A 



, A 



R2(t) = -p / Vn A -Vn-u/ I p A n A -An, 

Jo JT n JO JT n 

R*(t) = v / / p A n A V • (Vn Vn) + v / / (w • V)n A • An 

Jo JT n JO il" 

-v l I {u X ■ V)n A ■ An-v l [ (u ■ V)n • An A , 

JO JT™ JO JT™ 

i? A (i) = i/ /" [ (u x ■ V)n A • n + v [ [ (u ■ V) 

70 if" JO JT™ 



n • n 



R x {t) = -uel f (An + |Vn| 2 n) • n A - vO f f (An A + |Vn A | 2 n A ) • n, 

JO JT™ JO JT™ 

R x (t) = -v6l f (An+ |Vn| 2 n) • |Vn A | 2 n A - vO f [ (An A + |Vn A | 2 n A ) • |Vn[ 2 n. 

JO JT™ Jo JT™ 



Now we estimate each term above. Firstly, we have from (|2,9|) that 



\R x (t)\ 



< u 



p x \u x \ 2 



2\ 2 



+ IMI 



P X \u X o\ 2 



Vo A - 1 



1 

2\ 2 



< IM 



p X \u x \ 2 



T" 



p X - 1 



2\ 2 



< c 



1 



V + i 

(9 A \JT™ / 

oo 

< c\-\ 

Then we estimate the second term R x (t) and get 



+ 



1 






( 


+ 1 


oo 



7 ip a -ii ; 

^JT™ 



IKIbf/ P x \u x \ 

\JT™ 

/ |Po A "l| 

JT™ 



T" 



p X - 1 



Vo A + 1 



(5.12) 



\R$(t)\ = 



JT™ 



77 < 

JO JT™ 



(\/p^n A - u) ® {^/p^u X - u) 

2 



Vu+ / / (p A - • [(« • V)u] 



< CX^ + C 

< cx- 1 + c I f 

Jo Jt 



JT" 
t 



p x u x - u\ 2 , 



p x u x - u) ■ V 



JT« 



|lt|" 



where we have used the incompressibility V • u = to give the following fact: 



t r / u,i2 




/o Jt™ V ^ 



/7 ^(^-V-vf^) + ft A A -v(^f 

JO JT" \ z J JO JT™ V ^ 

^(x/p 1 - l)u x ■ V 



JT™ 

' /■ \ I / /■' 

A 1 |2 



* f A. MP 
2 / ./o jt™ 2 



Ai2 



< c / / + c / / |p A | 

\J0 JT" / \J0 JT™ 



Next we estimate the third and forth term and get 

\R x s(t)\ = I f [ p"p\ 

JO jt™ 



< < mi/7 (77 ^ - rA " 

'0 JT" / \J0 JT" 



|^(t)| < m / / I(i-pV-a«| 

JO JT" 

< c/77 ii-p a i 2X : ' ' 

'0 JT™ 




JT™ 



Now we turn to estimate R x (t) as 



Rg(t) = v\ / p A (u A - V)n- An + zv / / p u • \7 



JT™ 



/0 JT" 

4n 



(5.13) 



(5.14) 



(5.15) 



u x \ 2 \Au\ 2 ) <CX-\ (5.16) 



■>A,,\ .. / / „A/..A r-w. \.. , .. / / „A,.A y- / I x7n l" \ / / ^ A a ,„A 



+ i/ / (u- V)n A • An ; 



jt™ 



— V 



n(u x ■ V)n A -An-uf f (u • V)n • An A 

/ f (p x - l)(u x • V)n • An + 1/ / / (u A • V)n • An + 1/ / / (p A - l) t — 
'o JT n JO Jf n Jo Ji n 2 

t r ft r ft f 



Vnl 



n(n • V)n • An ; 
_' n 



+^ / / (n • V)n A • An A -u {u ■ V)n A • An - v 

Jo Jj n Jo Jt 

u f I (p x - l)(n A • V)n • An + v f I (p x - l) t ^f- 
Jo JT n Jo JT n z 

-v I [ [(u x - u) ■ V](n A - n) ■ An + v [ [ [u ■ V)(n A - n)(An A - An) 
Jo Jt u Jo Jt™ 

v f l {p X -l){u X -V)n-/\n+v f [ p x ^f- 
Jo JT n JO JT n z 

[(\fp x ~u x - u) ■ V] (n A - n) • An - / / [(1 - s/p*)u x ■ V] (n A - n) • An 
io JT n Jo Jr n 

[ [ (Vn-V)(n A -n)-(Vn A -Vn), (5.17) 
lo Jj n 




where we have used the incompressibility V • u = in the last inequality. Then by the 
Cauchy inequality and the uniform estimates (12. 9p . we get 



|^(/)j < C[f I |p A -l| 2 Vf/7 |n A | 2 |Vn| 2 |An' 2X " : 
lo Ji n / \Jo Jt™ 



t r \h / ft r | Vn | 4 \ 5 



+c( I I 



'O JT" / \J0 JT n ^ 

l i 

+C||An|| 00 f /" / \^u x - u\ 2 \ ( f I |Vn A -Vn| 2 V 



x I / ft r _ \ I 



+||u A || oo ||An|| 0O / / |Vn A -Vn| 2 / / |1 - vV 
+C||V«||oo / [ |Vn A - Vn| 2 

Jo JT™ 

< C{\- 1 + \- 2 ) + C ( [ \^fp x ~u x -u\ 2 + C (\\Vn\\l + l) I I |Vn A -Vn| 2 . 
Jo Jj n Jo Jf n 

(5.18) 

Noting that |n| = 1 and |n A | = 1 in Qt, we have the following two facts: 

V(|n A | 2 ) = V(|n| 2 ) = 0, (5.19) 

and 

(An + |Vn| 2 n) • n = (An A + |Vn A | 2 n A ) • n A = 0. (5.20) 
With the help of (|5.19p . we have 



R*(t) = v ( n ^ . v)n A • (n - n A ) + z> / / (u ■ V)n • (n A - n) 

Jo Jf n Jo Jj n 
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/' / / [(u x - u) ■ V]n A • (n - n A ) + v \ \ (u ■ V)(n - n A ) • (n x - n), 

(5.21) 



then following by the Cauchy inequality, we obtain 



\Rl{t)\ < C\\Vn> 



t r \ 2 / rt 

\u x -u\ 2 




JT n 




\n x — n\ 2 



o JT n 



I ( i 
\n x -n\ 2 ) 2 ( I I |Vn A - Vn| 2 



< C 



JT™ 



u x -u\ 2 + C 



rt 

/0 JT n 




JT" 



|n A -n| 2 + C 




JJ n 



|Vn A -Vn| 2 . (5.22) 



+C||«||oo 

/7 

JO JT n 

According to (I5.20p . we have 

R X (t) = -vQ [ [ (An + |Vn| 2 n) ■ (n x - n) - v6 [ [ (An A + |Vn A | 2 n A ) • (n - n x ) 
Jo JT n Jo Jf n 



-vB [ [ [(An + |Vn| 2 n)-(An A + |Vn A |V)] • (r 
Jo Jt" 



n — n), 



(5.23) 



then we use the Cauchy inequality to get 



|-R A (t)| < ^ / / |(An+ |Vn[ 2 n) - (An A + \Vn x \ 2 n x )\ 2 + C I I \» - n 

4 Jq Jjn Jq Jjn 




A|2 



(5.24) 



Similarly, with the help of (|5.20p . we derive that 
R x (t) = 



vO f* [ [An+\Vn\ 2 n) -\Vn x \ 2 n x -v6 [ [ (An x + \Vn x \ 2 n x ) ■ \Vn\ 2 n 
Jo Jf n Jo Jf n 

vQ f I (An+ |Vn| 2 n) • |Vn A | 2 (n A - n) 
Jo Jj n 

u9 [ [ (An A + |Vn A | 2 n A ) ■ \Vn\ 2 (n - n x ) 
Jo Jt" 

t 



vQ I I (An A + |Vn A |^n A ) - (An + \Vn\ z n) -|Vn| 2 (n A -n) 
Jo Jt* 1 L J 

vQ \ / [(Vn A - Vn) • (Vn A + Vn)](An+ |Vn| 2 n) ■ (n A — n) 
Jo Jj n 



By the Cauchy inequality, we have 

|(An A + \Vn x \ 2 n x ) - (An + \Vn\ 2 n)\ 2 + C\\Vn\ 



4 .10 JT n 





(5.25) 



\n x - n\ 2 



JT" 



< 



+C||An+|Vn| 2 n|| 2 J||Vn A || 2 + ||Vn|| 2 / / |Vn A - Vn| 5 

+C f f \n x -n\ 2 
Jo JT™ 



\(An x + |VnYn A ) - (An + |Vn|^n 



,A|2„A 



2„\|2 



JT™ 



4'? 



+C||Vn||| / f \Vn-Vn x \ 2 + cf f 
Jo Ji n Jo Jt t 



\n x - n\ 2 . 



(5.26) 



u6 



In conclusion, combining all the estimates about R x (t) for i = 1, 2, . . . , 8, we have 

/ [\\v r p x ~ uX - u ? + 7tI™ a - n \ 2 + 7Tl VrjA - Vn l 2 + n A (>, 

+/i / / |Vn A - Vu| 2 + (k + u) / / |V-u A | 2 
jo ip Jo ip 

|(An A + |Vn A | 2 n A ) - (An + |Vn| 2 n)| 2 




JT n 



< 



Po u o - u o 



+ -|Vn A - Vn | 2 + n A (x,0) + CA" 1 



+C(||Vn||| + l) / f (|VpV -u\ 2 + |n A - n| 2 + |Vn A - Vn| 2 ) . (5.27) 



Noting the following facts that 

2 



' Po n o - ""o 



/ Jiy -u x +[ \u x - Uo \ 2 <c(\- 4 +\- 

Jjn v Jfri 



(5.28) 



and 



|u A -d 2 < 



yj~f?"u x 


-u X \ 2 + [ 




Jv* 




V / P A 'U A - it 2 







p A u A - u| 2 



(5.29) 

By (l5T27l) - (j5T29]l . the Gronwall's inequality and the fact that J Q ||Vn||| < C, we have 

||u A - u\\ 2 + ^||Vn A - Vn|| 2 < CA" 1 , (5.30) 



and then integrating (|5.27p over [0, t] to get 

I (ji\\Vu x -Vu\\ 2 + v9\\{An x + \Vn x \ 2 n x ) - (An+\Vn\ 2 n)\\ 2 ) <C\-\ (5.31) 
Jo 

Finally, we have from (jl.38|) c. and (|1.40p ^ that 

(n A - n) t - 6A(n x - n) = - [(u x - u) ■ V] n x - (u ■ V)(n A - n) 

+9[(Vn x - Vrt) : (Vn A + Vn)] • n x + 6\Vn\ 2 (n x - n), (5.32) 

then we have by the parabolic theory that 

\\n x -n\\ 2 < C||Vn A ||^||u A - n|| 2 + 011^11^1^^- Vn|| 2 

+C (||Vn A |& + \\VnfJ) ||Vn A - Vn|| 2 + C||Vn||^||n A - n|| 2 
< C\-\ (5.33) 
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Applying V to ([5.32p . we have 

V(n A - n) t - 6»AV(n A - n) 
= -[V(u x - u) • V]n A - [(u x - u) • V]Vn A - (Vu • V)(n A - n) - (u ■ V)V(n A - n) 
+fl[(VV - V 2 n) : (Vn A + Vn)] • n x + 9[{Vn x - Vn) : (V 2 n A + V 2 n)] • n x 
+6»[(Vn A - Vn) : (Vn A + Vn)] • Vn A + 20(Vn : V 2 n) • (n A - n) + 6»|Vn[ 2 (Vn A - Vn), 



JO 

Now we complete the proof of Theorem 12.31 □ 
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(5.34) 



and then (|fT30j) . (|fT3T|) and (1031) imply 




(5.35) 
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